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Abstract 

We formulate a general framework for describing the electromagnetic 
properties of spacetime. These properties are encoded in the 'constitutive 
tensor of the vacuum', a quantity analogous to that used in the description 
of material media. We give a generally covanant derivation of the Fresnel 
equation describing the local properties of the propagation of electromag- 
netic waves for the case of the most general possible linear constitutive 
tensor. We also study the particular case in which a light cone structure is 
induced and the circumstances under which such a structure emerges. In 
particular, we will study the relationship between the dual operators de- 
fined by the constitutive tensor under certain conditions and the existence 
of a conformal metric. Closure and symmetry of the constitutive tensor 
will be found as conditions which ensure the existence of a conformal met- 
ric. We will also see how the metric components can be explicitly deduced 
from the constitutive tensor if these two conditions are met. Finally, we 
will apply the same method to explore the consequences of relaxing the 
condition of symmetry and how this affects the emergence of the light 
cone. 
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1 Introduction 



1.1 Motivation 

In Einstein's theory of gravity, General Relativity (GR), the fundamental vari- 
able describing the gravitational field is the spacetime metric. Einstein's equa- 
tions are partial differential equations which, together with suitable boundary 
conditions, determine the metric field from the energy-momentum distribution 
of matter. 

In GR, the metric tensor fixes all properties of spacetime: its metric prop- 
erties (allows to define lengths), its affinc structure (allows to compare object 
at different points, by parallel transport) and also its electromagnetic properties 
(specifies how electromagnetic waves propagate and, in particular, defines the 
light cone at each event). 

In searching for a better understanding of the structures which can be used 
to describe spacetime it is interesting and useful to distinguish each different role 
that the metric plays in GR. This approach has been followed before regarding 
the affine properties of spacetime. As a result, it is now clear that these proper- 
ties are in general described by an affine connection, which is an object logically 
independent of the metric. It has been useful to consider the metric and the 
affine connection as independent fields describing different aspects of spacetime. 
In this way, it was possible to formulate more general, so-called metric-affinc 
theories of gravity and spacetime which have the potential to describe physics 
beyond GR. In this more general framework, Einstein's theory is recovered as 
a particular case when certain additional conditions (vanishing of torsion and 
non-metricity) are satisfied. When this happens the affine connection turns out 
to be a function of only the metric tensor, as in GR. On the other hand, the 
more general formalism allows us to search additionally for alternative points 
of view as, for example, the possibility to formulate a purely affine theory of 
gravity, in which the affinc connection is considered to be the fundamental field 
and the metric would be a derived quantity. 

One of the aims of the present work is to extend this line of thought also 
to the description of the electromagnetic properties of spacetime. For this, it is 
necessary to recognize as clearly as possible the role which the metric and affinc 
structures play for the formulation of classical electrodynamics. It will turn out 
that the basic structure of electromagnetic theory, when properly formulated, 
is independent of metric and affine concepts. The electromagnetic properties of 
spacetime are encoded in the so-called constitutive tensor of spacetime, which 
is a quantity analogous to that used for the description of electrodynamics in 
material media. We consider the constitutive tensor of spacetime as indepen- 
dent of metric and affine structures, in the same way as the affine connection 
was considered to be independent of the metric, and develop the formalism for 
the description of the electromagnetic properties of spacetime. Not surprisingly, 
structures richer than in GR are found. In particular, we study some aspects of 
wave propagation in the case of an arbitrary constitutive tensor and derive the 
Fresnel equation describing the local properties of the propagation of light. In 
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general, instead of a light cone structure related to quadratic forms, a general 
quartic structure is found. The particular case of Einstein's theory corresponds 
to a particular choice of constitutive tensor, completely determined by the (con- 
formal) metric. 

We study the conditions which must be satisfied such that the electromag- 
netic properties of spacetime are compatible with a light cone structure. Finally, 
it is also possible (similarly to the case of purely affine theories) to consider the 
constitutive tensor as a more fundamental field describing spacetime and the 
conformal metric as a derived quantity. 

That the conformal metric can indeed be considered as a consequence of 
the electromagnetic properties of spacetime should be intuitively clear. The 
causal properties of spacetime, described by the conformal metric in GR, are 
of an electromagnetic nature. Light rays are the physical realization of null 
vectors fields, null hypcrsurfaces are the regions through which electromagnetic 
waves can propagate. From an operational point of view it is clear that most of 
our knowledge of spacetime is extracted from the properties of electromagnetic 
fields and how they propagate. Today, we send electromagnetic waves towards 
(reflecting) objects to measure distances. We measure velocities by means of 
Doppler-shifted electromagnetic waves, etc. The properties of spacetime itself 
are then defined as those independent of the specific configuration of the fields 
and interpreted as describing the underlying 'substrate' on which we make mea- 
surements. 

The idea of considering the spacetime metric as a secondary, derived field, is 
not new. Already in 1921 Eddington studied a 'purely affine gravity' model, in 
which the metric is defined in terms of a symmetric affine connection, see |64j . 
Thus, in Eddingston's theory, the causal properties of spacetime are derived 
from its affine properties. Other models have proposed to replace the metric by 
a trio of self-dual 2-forms [3 [9] . In [10] the metric is obtained from a solution 
of a theory formulated only in terms of a SL(2,C) connection, a tetrad, and a 
scalar density. More recently, it has been shown by Barcelo et al. [H [3] that 
an effective metric can be derived (defined) for almost any lagrangian theory 
of scalar fields, provided one considers perturbations of the fields around some 
background configuration. If the theory depends of a single scalar field, then 
a metric can be uniquely defined (so that the equation for the field perturba- 
tion can be written as a Klein-Gordon equation with respect to that metric) 
[2 . If more scalar fields are involved, multiple metrics can be introduced in 
general, which can be described as refringence (birefringence, trirefringence, ...) 
in the sense that different fields of the theory would 'see' different effective met- 
ric structures [3]. In these models, the causal properties of space are again a 
manifestation of the dynamics of some more fundamental fields. 

In a related context, there has been renewed interest in recent years for so- 
called 'analogue models of gravity', see [78]. These models are based on the 
results that some condensed matter systems, as, for instance, a moving fluid 
with acoustic perturbations, a moving dielectric with light, or a moving super- 
fluid with quasiparticlcs, can be described in term of some 'effective metric'. 
The metric is also here a derived, secondary object which turns out to be useful 
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for a geometrical description (a la GR) of the system properties and depends 
on the more fundamental degrees of freedom of the system. The interest in 
these analogue models lies in the possibility to construct systems in which kine- 
matical properties of physics in curved space could be simulated and tested in 
the laboratory. In [55], for instance, the possibility of constructing a 'dielectric 
black hole analogs' is discussed, i.e., dielectric materials in which the effective 
metric describes the analog of an event horizon, see also |40[ 177] . Dynamics, on 
the other hand, is not likely to be simulated, since in general the dynamics of 
the effective metric can be completely different to that imposed by Einstein's 
equations in GR. 

My goal here is not to propose a concrete alternative spacetime theory, 
but the more modest one of studying a more general framework describing 
classical electrodynamics in which a metric tensor is not assumed as a basic 
fundamental field from the very beginning. This will raise the question on which 
fundamental structures classical electrodynamics can be based without using a 
metric. The framework will be called pre-metric (or sometimes metric-free) 
electrodynamics, and it shares many features and analogies with the theory of 
electrodynamics in a material medium. By studying such a general framework, 
the conformal properties of spacetime, i.e., the conformal spacetime metric with 
its nine independent components, can eventually emerge as an special case under 
certain particular circumstances. The study of the conditions under which a 
conformal metric structure is induced is also one of the subjects of this work. 
This framework could then be useful to develop a generalized theory of spacetime 
in which the properties of vacuum are treated in analogy to a material medium 
and in which GR with its causal structure could be recovered in suitable limiting 
cases. 

The axiomatic approach to classical electrodynamics presented here has been 
developed from the original ideas of Kottler [36] and van Dantzig [76]. They 
seem to be the firsts who recognized that the fundamental structure of Maxwell's 
equations is independent of the metric and affine structure of spacetime. In what 
concerns the derivation of the spacetime metric from linear electrodynamics, in 
which we are interested here, Peres [5B], already in 1962, wrote: 'It is therefore 
suggested to consider the electromagnetic field as fundamental, and the met- 
ric field only as a subsidiary quantity'. This same idea was also developed by 
Toupin [71] and by Schbnberg |62j who showed that a conformal metric structure 
is induced assuming that the 'constitutive tensor' defining the 'spacetime rela- 
tion' between electromagnetic excitations and field strengths (see section 12. 6|) 
satisfies the conditions of symmetry and closure (section 0} . Jadczyk [3T] also 
showed that a spacetime metric can be introduced under the above mentioned 
conditions. However, only very recently and explicit derivation of the induced 
conformal metric has been given by Obukhov and Hehl [52"] . 

We formulate a general framework for electrodynamics in an arbitrary lin- 
ear electromagnetic medium^. These developments are useful for three different 

1 We use the word 'medium' in a general sense in order to refer to any 'arena' on which 
electromagnetic phenomena could take place. In particular, we include the vacuum as a 
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aspects at least. First, they make the fundamental structures of classical elec- 
trodynamics more transparent. Second, they may provide a basis for a deeper 
understanding and for generalized models of spacetime and its electromagnetic 
properties. This includes, for instance, the application of the formalism to study 
test theories, see |21) and references therein. Third, the formalism developed 
in sections [5] and [3] can also be interpreted and used in optics as a general co- 
variant theory of electrodynamics in inhomogeneous, anisotropic, and in general 
dissipative material media. 

Finally, we study the particular case in which a light cone structure is induced 
and the circumstances under which such structure emerges. In particular, we 
will study the relationship between the dual operators defined by the constitutive 
tensor under certain conditions and the existence of a conformal metric. Closure 
and symmetry of the constitutive tensor will be found to be the conditions which 
ensure the existence of a conformal metric. We will also see how the metric 
components can be explicitly derived from the constitutive tensor if these two 
conditions are satisfied. We will also give an alternative, simpler, and more 
physical derivation of the metric based on direct use of our general results about 
the Fresnel equation describing the local properties of light propagation. Finally, 
we will apply the same method to explore the consequences of relaxing the 
condition of symmetry and how this affects the emergence of the light cone. 

1.2 Maxwell-Lorentz equations in 3-vector notation 

We start with electrodynamics within the framework of Special Relativity (SR), 
but in 3D formalism. 

In the usual 3-dimensional vector notation, V := (V x , V y , V z ), the Maxwell- 
Lorentz equations in integral form read: 

f D • dS = [ pdV, <f H • dr = [ J • dS + — ( [ D • dS ) , (1) 

JdV JV JdS JS ™ \Js J 

L Bds =°' L'-"—W. B -")- (2) 

Here D is the electric excitation (historically called 'electric displacement'), p 
the electric charge density, H the magnetic excitation (historically called 'mag- 
netic field'), J the electric current density, B the magnetic field strength and 
E the electric field strength. The integrals above are defined over arbitrary 
volumes V with boundary dV and over arbitrary surfaces S with boundary 
dS, respectively. The corresponding volume elements of the three-, two-, and 
one-dimensional regions are denoted by dV, dS, and dr. Finally, the dot ■ de- 
notes the 3-dimensional scalar product of vectors, which is a metric-dependent 
operation. 

particular electromagnetic medium. On the other hand, 'material' or 'material medium' will 
be used to refer to media with a known (atomic) substructure as, for instance, crystal, liquids, 
etc. 
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Equation (jTJi) summarizes the Gauss law. It implies that the field lines 
defined by D can be open, the ends of which are located at points where charge 
is located. Equation $2&) is usually interpreted as expressing the absence of 
magnetic monopoles in nature0. It implies that the magnetic lines defined by 
B must be closed. Equation ^p) summarizes Faraday's induction law ('a time- 
variation of an magnetic field induces an electric field'). 

The field strengths E and B are operationally defined by means of the 
Lorentz force law. On test charges, the force density is given by 

F = pE + JxB. (3) 

Here x is the 3-dimensional vector product, also a metric-dependent operation. 
The Maxwell equations are completed by the relations 

D = e E, H = — B. (4) 
Mo 

The constants Eq and /xn are called permittivity and permeability of vacuum, 
respectively. The speed of light, i.e., the speed with which electromagnetic 
perturbations propagate, is given by c = -j==. 

As it is well known, the Maxwell equations ([T]) and (J2j) can be written in 
differential form as 

V-D = p, VxH = J + — , (5) 

V-B = 0, VxE=- — . (6) 

The inhomogeneous Maxwell equations ((SJ) are such that the conservation of 
electric charge is automatically satisfied, i.e. 

f + V ' J = °- ^ 



Actually, Maxwell completed the electromagnetic equations known at his 

3D 

m 



time by adding the 'electric displacement' term such that the resulting equa- 



tions were consistent with charge conservation. 



1.3 Electrodynamics in a material medium 

It is well known that a macroscopic description of electromagnetic phenomena 
inside a material medium (treated as a continuum) can be achieved using macro- 
scopic Maxwell equations which are of the same form as (J5J) and ([6]), but were 
now 

• D and H denote macroscopic field strengths related to the sources p and 
J, which are now external charge and current densities, respectively, and 

2 For recent (unsuccessful) searches for magnetic monopoles, see |22l HI |32] and references 
therein. 
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• the relation between (D, H) and (E, B) is now not given by ((31) but by a 

constitutive relations 

D = D [E, B] , H = H [E, B] , (8) 

which contain the information of the particular electromagnetic properties 
of the medium under consideration. 

Among the many possible particular cases (nonlocal constitutive laws, ...) 
we recall here the case of a linear anisotropic medium for which the constitutive 
relations are usually written, in components (see for instance |38| . page 313), as 

D a = £ab Eb + a a bBb, H a = fl ab Bb + PabEb- (9) 

Here a, 6, . . . = 1,2, 3, e a b is the permittivity tensor, and [i a b the permeability 
tensor of the medium. The tensors a and (3 describe magneto-electrical prop- 
erties. For examples of and further discussions on magneto-electric media, see 

Consider the case in which a a b and (3 a b vanish. A non-magnetic medium 
corresponds to the subcase in which ii a b = fJ>o5 a b so that the vacuum relation 
IJDd) holds. Interesting properties of non-magnetic anisotropic media regarding 
propagation of plane electromagnetic waves include [38] : 

• In general, for a triaxial crystal (i.e., when the three eigenvalues of e a b 
are different), a fourth order Fresnel equation determines the dispersion 
relation of plane waves. 

• The wave vector k and the ray vector s are in general not parallel: k is the 
vector normal to the wave front, s is the direction of energy propagation. 

• Birefringence in uniaxial crystals (two eigenvalues of e a b are equal): the 
Fresnel equation factorizes into two quadratic factors, one isotropic factor 
corresponding to 'ordinary' waves and an anisotropic one corresponding 
to 'extraordinary' waves. 

Effects analogous to the above mentioned will be discussed in section [3J 
Other examples of local constitutive laws include: 

• Double refraction induced by an electric field. This so-called Kerr effect 
can be induced if an isotropic material is placed in a constant electric field. 
The electric field breaks the isotropy of the medium producing a change in 
the dielectric constant, leading to effects similar to those observed in uni- 
axial crystals. This effect can be described by the nonlinear constitutive 
law corresponding to the following dielectric tensor: 

£ab = e 5 ab + aE a E b , (10) 

with some constant a. 

• Magneto-optical effects: The dielectric constant depends on the magnetic 
field strength H. 
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vector/scalar 


P 


J 


H 


D 


E 


B 


p-form 


P 


3 


H 


V 


E 


B 


P 


3 


2 


1 


2 


1 


2 



Table 1: Correspondence between vectors and exterior forms. 



1.4 Maxwell-Lorentz equations in 3+1 exterior form no- 
tation 

From the Maxwell-Lorentz equations in their integral form, see (JJ) and ([2]), 
one can see that the different fields appear associated with integrals over re- 
gions of different dimensionality. For the sources, we see that p is integrated 
over 3-dimensional regions (volumes) while J is integrated over surfaces, i.e., 
2-dimensional regions. The excitations H and D are integrated over 1- and 
2-dimensional domains, respectively. Finally, the field strengths E and B ap- 
pears in the Maxwell-Lorentz equations under 1- and 2-dimensional integrals, 
respectively. From the theory of exterior forms, see [16] , we know that a p-form 
is the natural object to be integrated over a p-dimensional domain. This means 
that the Maxwell-Lorentz equations can be reformulated in terms of exterior 
forms according to the identifications of table [1] In terms of exterior forms, the 
Maxwell-Lorentz equations ([I]) and ([2]) are naturally expressed as 

dV = p, dH=^-V + j, (11) 
d 

dB = 0, dE+—B = 0. (12) 

Here d = dx a A d a denotes the 3-dimensional exterior derivative (we will use d 
for the 4-dimensional exterior derivative) . Note the internal consistency of these 
equations, since d increases by one the rank of the exterior form on which it is 
applied. Conversely, taking the equations (fTTjl and (fT2|) as starting point, one 
can derive the corresponding field equations in terms of field components. We 
decompose each form as follows: 

P = 3~! Pabc ^ A ^ A ^ = A ^ A 
3 = \ 3ab dx a A dx b = ^ e abc f dx a A dx\ (14) 

n = n a a- v = \ v ab dx « a ^ = \ ^ a dx\ (15) 

E = E a dx a , B = ^B ab dx a Adx b = ^e abc B c dx a A dx b , (16) 

Here e a b c = £[ a bc], with £123 = 1, is the 3-dimensional Levi-Civita symbol. Then, 
substituting (p~5]) — i|T^|) into the equations (TTTj) and l|12p . one directly finds, see 



10 



d a v a = P , e abc d b n c = JV + f, (17) 

d a B a = 0, e abc d b E c + ^B a = 0, (18) 

which generalize §5§ and ©. Note that T> a and H a are 3-dimensional vector 
densities. The formulation in terms of exterior forms has the advantage that the 
corresponding p-forms are independent of the 3-dimensional coordinate system 
used. In other words (fTTj) and (fT2"|) , and therefore also (fTT)) and (fT5)) . are valid 
not only in cartesian coordinates but in any 3-dimensional coordinate system. 
Actually, if one considers the exterior forms in (jT3])— (jTHJ) as basic field variables, 
as we do here, then equations (|11[) and (|12[) are independent of any metric or 
affine structure of the 3-dimensional space. 

1.5 Poincare covariant Maxwell equations 

As usually shown in textbooks on SR, the Maxwell equations (0 and (J6j) can be 
written covariantly under Poincare transformations by defining the field strength 
Fij = —Fji and the excitation Hij = —Hji = t, x, y, z) by 

(19) 

(20) 
(21) 
(22) 

Additionally, the electric current 4-vector density J 1 is defined as 

J l :=p, J X :=J X , J y :=J y , J z := J z . (23) 
Then, we can write ([5]) and (J6]) as 

^ kl d : H kl = J\ JWdjFu = 0, (24) 

where e Ijfez is the Levi-Civita symbol. The relation between Fij and Hij , namely 
equation (j4|), is translated into 

H v = Ix^iJkl V km r] ln F mn , (25) 
^ V / i0 

where ifi are the components of the Minkowski metric in cartesian coordinates: 
r]ij = diag(c 2 , — 1, — 1, — 1), rfi = diag(c -2 , — 1, — 1, — 1). Inserting (|2"B")) into 
(|24|) we obtain (always in cartesian coordinates) the inhomogeneous Maxwell 
equations in a form which is often used in SR, namely 

x [^djF ij =J\ (26) 
V Mo 



Ftx '■— 


-E x , 


F ty :■ 


= -E\ 




:= -E z 


TP 

r X y 


:=B\ 


TP 

r yz 


:= B x , 


F ■ 

± zx • 


= B v , 


Htx '■ 


= H x 1 


H t y 


:= H y , 




■=H\ 


H xy 




Hy Z 


:= D'\ 


H zx 


:= D y . 



11 



1.6 Poincare group invariance and 'natural' invariance 

It is well known that the Maxwell equations (|24[) and (f2"5|) are form invariant 
under Poincare transformations of the form 

x l -> x 1 ' := K l \x l + a\ (27) 

with A t i 6 50(1,3), i.e. satisfying 

A<A^' =77 <Y , a* eft. (28) 

This form invariance of the Maxwell equations means that the physics of the 
electromagnetic fields remains the same on rotated frames and in frames moving 
with constant velocity with respect to each other (i.e. boosted frames), as well 
as at different points in spacetime. 

However, a more careful analysis shows that the field equations (|24|) are ac- 
tually form invariant under any coordinate transformation x l — > x % = x % (x 1 ), 
provided Fij, and Hij are considered as components of a second order (anti- 
symmetric) tensor field, and J 1 as the components of a vector density field of 
weight +1. This feature of (|24[) is sometimes called 'natural invariance' of the 
Maxwell equations, see |59j . chapter 3, for an extended discussion. The natural 
invariance shows, on the other hand, that the physical information about the 
equivalence of inertial frames under Lorentz transformations is contained exclu- 
sively in the spacetime relation between excitation and field strength (123)) . This 
is a property of the vacuum. Furthermore, we note that the Minkowski metric 
only appears in (f25|) . Looking at (|25|) it is then clear that this invariance of the 
vacuum is a direct consequence of the form invariance of the Minkowski metric 
with respect to Poincare transformations, see (f28|) . In section 12. 7. II we will give 
a more general definition of symmetries of a given electromagnetic medium. 

1.7 Maxwell equations in curved spacetime 

The traditional recipe to construct a theory including the interaction with the 
gravitational field is to 'replace partial derivatives by covariant derivatives', also 
called 'minimal coupling' to the gravitational field. This procedure ensures that 
the resulting equations are covariant under an arbitrary change of coordinates, 
see, for instance, [13] for more details and examples. The procedure is not free 
of ambiguities when one considers the electromagnetic potential as fundamental 
variable, since then the Maxwell equations are of second order and a 'normal or- 
dering' problem appears when applying the recipe above. This is a consequence 
of the fact that covariant derivatives do not commute in a curved spacetime. 

A not widely recognized fact is, however, that in the case of electrodynamics 
the mentioned recipe is completely unnecessary since the Maxwell equations, 
when properly formulated, are 'naturally covariant', as we have seen in section 
11.61 For the (rather trivial) transition from a Minkowski space to a Riemannian 
space with metric g one just needs to replace the spacetime relation by 

Ha = <[^e m V\g~\9 km 9 ln F mn , (29) 
V Mo 
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with g := det {gij), so that F and H are tensors, as required by the natural 
invariance of (|24|) . At every event one can find Riemann normal coordinates 
which will reduce the metric to its minkowskian form, i.e. gij = rjij and then 
(f!?5|) reduces to the form in agreement with the equivalence principle. 

Equation (|29[) can be written in terms of the Hodge dual operator * of the 
metric g, namely 



The Maxwell equations in a Riemannian space are therefore, in terms of F 
and g: 



The fact that the Maxwell equations ([M)) are naturally invariant and that the 
metric structure of space does not enter in their formulation is not an accident. 
It is a consequence of the fundamental property that the basic structure of 
electrodynamics can be derived from counting procedures of charge and magnetic 
flux, with do not require a metric (nor an affine) structure of spacetime. This 
properties will be further clarified in the axiomatic approach of section [5J 

2 Electrodynamics on an arbitrary 4D-manifold 

In this section, we would like to present an axiomatic construction of classical 
electrodynamics which intends to be as general as possible. Structures are only 
introduced when they are indispensable for the development and not earlier 
than necessary. This approach will provide us a very general framework which 
can then be applied to many different particular cases. 

We model spacetime as a smooth 4-dimensional manifold X which, at least 
in some neighborhood, admits a foliation into 3-dimensional submanifolds, pa- 
rameterized by a monotonic 'time' parameter a, see figure [TJ 

2.1 Charge conservation 

Probably the most important and defining property of electromagnetic theory 
is the experimentally well tested fact of charge conservatiord. A basic property 
of electric charge is that it is an additive quantity which can be distributed in 
space, i.e., it is an extensive quantity. In nature, electric charge is know to be 
quantized, its fundamental quanta, the electric charges of the quarks are ±e/3 
and ±2e/3, and those of the leptons ±e, where e denotes the elementary charge. 
In a classical field theory, we describe the distribution of an extensive quantity 
in terms of a current density, containing the information of how many charges 
are distributed in spacetime, and how they move. 

3 For experiments testing charge conservation/violation see, for instance, |49l 167] 




(30) 





(31) 
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Figure 1: Foliation of spacetime: Each hypersurface h a represents, at a 'time' 
<7, a 3-dimensional submanifold. 

Consider the total charge within some 3-dimensional region. This total 
charge should be independent of the orientation of any coordinate system used. 
Then, the electric charge distribution is described by the twisted electric current 
3-form J, see [16] for definitions. For a 3-dimensional region S3 of a hypersur- 
face h a one can interpret J s J as the total charge contained in S3. On the 
other hand, if the 3-dimensional region is of the type — S 2 x [ai, a 2] one can 
interpret J" n J as the total charge crossing the 2-dimensional surface S 2 during 
the 'time' interval [a\, 02]. Consequently, the 3-form J carries the dimension of 
charge, i.e. [J] = q. 

Note that no concept of distance or parallel displacement, i.e., no metric or 
connection are necessary to define the concepts of charge and charge current 3- 
form specifically. Of course, quantities like 'charge per unit volume' and 'charge 
per unit area and unit time 1 are useful after one provides prescriptions for what 
'unit volume', 'unit area' and 'unit time' are. However, the latter are concepts 
not needed to describe, e.g., how many electrons and therefore how much charge, 
are contained in a certain region. Clearly, the total charge is independent of the 
unit in which volume is measured. 

The components Ja/37 of the electric current 3-form with respect to some 
coframe basis i? a , a, (3, . . . — 0, 1, 2, 3, are determined by 

J=~ J a ^ a A & A (32) 
If one associates a dimension I, in the sense of a segment to the coframe #, i.e. 

4 i.e., a one dimensional extension on the manifold, not in the sense of a unit of length. 
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Figure 2: Local conservation of charge: Each worldline of a charged particle 
that enters the finite 4- volume Q± via its boundary dQ^ has also to leave f^4 . 

[d] = I, then the components J a f3~f carry the dimension [J Qi g 7 ] = ql~ 3 . 

The conservation of electric charge is then expressed as the vanishing of the 
integral 

I J = 0, V n 4 , (33) 

i.e., for any 3-dimensional boundary of a 4-dimensional region ft^. In particular, 
for a region — £3 x [o~\ , 02] , the integral conservation law (|33p requires the 
balance between the change of the charge in the region £3 during the interval 
[ci, 1^2] and the flux across its 2-dimensional boundary 8^3. 

Since the region ^4 in (|33|) is arbitrary, the Stokes theorem tells us that the 
current 3-form must be closed, i.e. 

dJ = 0. (34) 

Now, according to the de Rham theorem, see for instance [16] . the current 
3-form is not only closed but also exact, i.e., it can be derived from some 2-form 
H by exterior derivation, if all its integrals vanish over 3-dimensional regions 
without boundaries. Under these assumptions, this means that there must be 
a 2-form H such that 

dH = J. (35) 

The twisted 2-form H is called the electromagnetic excitation and carries dimen- 
sion of charge, [H] = q. As we will see, in a 3+1 decomposition its components 
can be identified with the metric-free generalization of the electric and magnetic 
excitations. 
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However, the conditions above are not enough to uniquely define the 2-form 
H, since a 'gauge' transformation 

H -> H' :=H + <M> (36) 

leaves (|35p invariant, for an arbitrary twisted 1-form <5 . 

Equation (j3"5)) arc thus the inhomogeneous Maxwell equations (4 equations). 
More than defining the values of H, the above arguments show that the inho- 
mogeneous Maxwell equations must be of the form (|3"5|) , since this is the only 
kind of field equation which are compatible with electric charge conservation. 

A single electromagnetic excitation will be picked out by the requirement 
that H = for F = for the spacetime/constitutive relation, see sections 12.51 
and [231 



2.2 Lorentz force 

We assume now that the concept of force density is known from mechanics and 
use it to define the electromagnetic field strength 2-form F, as usual, as force 
per unit charge. We define F by means of 

fa =: (e a \F) A J, (37) 

where e a is a framed and f a G A 4 (X) are the corresponding components of the 
force density covector-valued 4-form in that fram€0. From the definition ([37} F 
is a untwisted 2-form, i.e., an intensive quantity. It describes how strongly the 
electromagnetic field acts on test currents. The dimension of F is [F] = h/q, 
with h denoting the dimension of an action. Equivalently F has dimension of 
magnetic flux. The definition above is a very restrictive one, since it tells us 
that the force on test charges is determined by only the six independent com- 
ponents of the 2-form F, instead of the 16 independent quantities that a linear 
relation between force density and current density would admit in principl^]. 
This assumption is part of our axiomatics, and it is suggested by the fact that 
the excitation is described by a field with 6 independent components and one 
therefore expects the field strength to have the same number of independent 
components. We know, of course, that this choice is reasonable since we know 
that in Maxwell's theory, the Lorentz force is determined by the six independent 
components of the electric and magnetic field strengths. 

It is true that in order to have a complete predictive theory, one still has 
to specify the relation between velocities and momenta of test currents. This is 
required in order to be able to predict the evolution of test currents in a given 
electromagnetic field. This relation includes the metric tensor in the known case 
of GR. However, we are interested here in the structure of the general electro- 
magnetic theory, and not specifically in the mechanical 'constitutive relation'. 

5 This is again a metric-independent quantity. A frame is just a basis of the tangent space. 

6 Remember, in classical mechanics force f; = Q~ and momentum p; = are both 

ox 1 or 

covectors. 

7 Consider for instance a relation of the form f a = ty a A J. Then ty a is a covector valued 
1-form, and has therefore 16 independent components. 
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2.3 Magnetic flux conservation 



The next step in our axiomatic construction is to find conditions for the elec- 
tromagnetic held strength to satisfy. The natural operation that can be done 
with a 2-form is to integrate it on a given 2-dimensional region. If we assume, 
in analogy to conservation of charge, that 

I F = 0, (38) 

for an arbitrary 3-dimcnsional submanifold Q3, then Stokes' theorem provides 
us with a differential equation for F, namely that the held strength must be 
closed, 

dF = , (39) 

and (at least in some neighborhood of each event) exact, i.e. F — dA. The 
untwisted 1-form A is then the electromagnetic potential. 

We take (|38|) . or equivalently (139|) as third axiom. It represents the ho- 
mogeneous Maxwell equations (4 equations) and expresses the conservation of 
magnetic flux. In general, magnetic flux is not quantized, as it is the case of elec- 
tric charge. In type II superconductors, however, as, e.g., in Niobium, quantized 
magnetic flux lines are possible. 

This formulation does not admit magnetic charges (i.e. magnetic monopoles) 
in a natural way. Since the held strength F is by its very definition an intensive 
quantity, a hypothetical magnetic charge density p m , such that dF — p m , would 
necessarily also be an intensive quantity (an untwisted 3-form) , quite in contrast 
to the extensive nature of any charge-like quantity like, e.g., electric charge, 
energy-momentum, all twisted 3-forms. 

Locally at least, see above, the homogeneous Maxwell equation (l39|) imply 
that the field strength F can be derived from a untwisted 1-form A, the elec- 
tromagnetic potential, such that 

F = dA. (40) 

Of course, A is only determined up to a gauge transformation 

A—y A' = A + d&, (41) 

for any untwisted 0-form (i.e. scalar) "J. 

The electromagnetic potential carries the same physical dimension as the 
electromagnetic held strength, namely [A] = h/q. 

2.4 'Space-time' Decomposition 

Consider a foliation characterized by the monotonic 'time-like' parameter a. 
Consider also a vector field n not lying on a a — const, surface, i.e. n\da =/= 0, 
so that it can be used to 'evolve' the folia, see figure [1] In particular, one can 
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Figure 3: Conservation of magnetic flux in spacetime: For an arbitrary 3- 
dimensional integration domain f^, the integral <f F vanishes. 

rescale the vector field to make it fulfill the normalization condition n\da = 1, 
which turns out to be useful. 

We will call any p-form to '3-dimensional with respect to n' if n\cu = 0. The 
'time-like' components of a '3-dimensional form with respect to rC in a basis 
adapted to the folia (eo ~ n) vanish. Then, given any p-form 4* one can define 
the p-form and the (p — l)-form both of which are 3-dimensional with 
respect to n in the above sense, as 

* := n\ (da A 4") , (42) 

and 

4'_ L :=nj4'. (43) 

For p = 4 we have J> = 0. Similarly, for p = 0, 4 f j_ = 0. 

These two '3-dimensional' quantities $ and 4'j_ contain the complete infor- 
mation of the original form 4'. The latter can be written as 

4 1 = da A *j_ + (44) 

One may also call j£ and da A 4 f j_ the 'transverse' and 'longitudinal' parts 
of 4 1 with respect to n, respectively, see |27) . 

From their definition, 4_ and 4 f j_ satisfy the following properties 

• x of a product, 

(*A$)_l = * ± A$ + (-1) p *A$_l, (45) 
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• _ of a product, 

g A $ = g A (46) 

• Lie derivative £„ and ± commute, 

(£„*) ± = £„*l, (47) 

• Lie derivative £„ and _ commute, 

£njg = (48) 

Additionally, it is useful to define a 3-dimensional part of the exterior derivative, 
d, such that for any p-form ^> 

d^ :=n\ (da A*). (49) 
It has the following properties: 

• d is 3-dimensional with respect to n, 

^# = nJ(d#)=Q, (50) 

• _i_ of a derivative, 

(d*) ± = £ n £-d*±, (51) 

• _ of a derivative, 

d* = rf£. (52) 

By means of this general decomposition procedure, we decompose the Maxwell 
equations ([33)) and (f3"9"| in terms of '3-dimensional' quantities. First decompose 
the fields J, iJ, F and A and introduce the notation 

j:=-J±, p:=J, H:=H±, V := H, (53) 

E := —Fj_, B:=F, tp := -Aj_, A := A. (54) 
The 4-dimensional quantities can be reconstructed according to 

J = -j A da + p, H = -H Ada + D, (55) 

F = E Ada + B, A = -ipAda + A. (56) 
Now one can take the inhomogeneous Maxwell equation (|35|) and find 

{dH - J) ± = (dH) ± -J± = C n H -dHj_ - Jj_ = C n V - dH + j, (57) 

and 

dH_J_= dH- J = dH- J = dV - p. (58) 
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Analogously, from the homogeneous Maxwell equations (f3"9")) we obtain 

(dF) ± = C n F- dF ± = C n B + dE, (59) 

and 

dF = dF = dB. (60) 

Thus, we find the Maxwell equations in a 3+1 decomposed form to be 

dV = p, dH = £ n V+j, dB = 0, dE + C n B = 0. (61) 

We can also decompose the law of charge conservation (|3"4")l . It is important 
to note that no information is obtained from dJ_ since this quantity vanishes for 
any 3- form J, as can be seen from the definition (j42|) . We then compute 

(dJ) ± = £ n J-dJ ± = CnP + dj, (62) 

so that charge conservation means 

C n p + dj = 0. (63) 

Finally, we decompose equation (|40|) . We find 

{dA) ± = C n A-dA x =£ n A+d<p, (64) 

and 

dA= dA = dA, (65) 

so that 

E = -dip - L n A, B = dA. (66) 

The decomposed Maxwell equations (|6"Tj) naturally generalize equations and 
(|12p . The Lie derivative C n is the natural generalization of the time derivative 
since it measures the change of the (integral of the) corresponding field between 
different folia, according to the displacement induced by the vector field n, see 
figure |U 

2.5 Measuring H 

In this section, we discuss a general procedure for measuring the excitation H. 
This can be done by using an idealized object, namely an ideal conductor. This 
special material is assumed to have the following two properties: 

1. In an ideal conductor, all charges are located on its surface. In other 
words, the ideal conductor is such that inside it no free charges can be 
found. If ^3 is the 3-dimcnsional region describing the conductor, then 
the free charges are all located on 9^3. 

2. In the 'rest frame' of the conductor the electric excitation T> vanishes. If 
at some event inside the conductor a volume element is spanned by the 
vectors e a , a, b, . . . = 1,2,3 and eo = n is a vector pointing in the zeroth 
independent 'time' direction, then we assume H a b := eb\e a \H = 0, or 
equivalently H = da A H, with n\do~ = 1 and e a \do~ = 0, see (|55j). 
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A' 




Figure 4: Under the displacement x l — > x n = x % + ari l (x) (a being an infinitesi- 
mal constant parameter) the integral J Q , ^ of ap-form 'J over the corresponding 

p 

p-dimensional mapped region fl' p is given by J n , ^ = J n ^/ + a J n L n ^. 



With these assumptions one can measure the excitation by use of 'Max- 
wellian double plates'. Consider two (uncharged) parallel plates made of an 
ideal conductor and locate them at the point P where the excitation should be 
measured. The field strength (whatever value it may have) will induce, through 
the Lorentz force, surface charges in the conductor. Separate now the plates and 
measure the charge Q induced in one of its surfaces. One can then integrate the 
inhomogeneous Maxwell equation over a volume with one side in one conducting 
plate and the other between the plates, see figure [5] Due to property one, in 
the limit of vanishing 'thickness' of f^, see fig. [5J one finds 

f dH= [ H = (e 2 J ei j H) p ~ (e 2 J e x \ H) cond , (67) 
Jq 3 Jdn 3 

and J — Q, so that 

(iJ 12 ) P - (H 12 ) cond = Q. (68) 

The second term on the left hand side of vanishes because of property 2 of 
ideal conductors. Therefore, the induced charge Q determines the component 
H12 of the excitation. Similarly, by orienting the plates differently one can mea- 
sure, e.g., -f/13. Furthermore, by changing the state of motion of the conductor 
(i.e., different 4- velocities) one can measure components of H which are, in the 
notation we are using, of the form (eijeoj-ff) = -Hoi- 
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Figure 5: Maxwellian plates. Here e\ and e2 span the surface element parallel 
to the plates. 

2.6 Constitutive relations 

As it is clear from the last sections, the Maxwell equations in their form (I35| 
and (|39|) are valid for any electromagnetic medium (i.e. including the vacuum) . 
They describe the general features of electrodynamics which follow from charge 
and magnetic flux conservation. In this form, Maxwell equations are valid, for 
instance, in vacuum in Special Relativity, but also if gravitational effects are 
included in the context of General Relativity or alternative theories, as for in- 
stance in those formulated in a general metric-affine spacetime. They are also 
applicable to any material medium, if one interprets J and H as macroscopic 
quantities. In other words, (|3"5|) and (|3l))) are general structures of electromag- 
netic theory. The additional structure that really defines the particular physical 
properties of the system under consideration, is the subject of this section. 

We saw that the Maxwell equations |35|) and (|39|) amount to only 8 inde- 
pendent equations for the 12 independent fields, H and F. A further relation 
between the excitation and field strength, i.e., a connection between the field 
generated by the charges and the field acting on test currents is therefore nec- 
essary in order to make the theory complete. Such a relation is referred to as 
the constitutive relation when corresponding to some material medium, and as 
the spacetime relation when describing the vacuum properties. The formalism 
presented here can thus be used in both situations. In the first one, the cur- 
rent 3-form describes the so-called 'external' ('free') currents flowing through a 
material medium. In vacuum, J represents the total current density. In any 
case, the relation between H and F describes the properties of the arena in 
which the electromagnetic phenomena of interest take place. Clearly, however, 
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Figure 6: The spacetime/constitutive relation connecting excitation and field 
strength. 



the spacetime relation is of more fundamental importance, since it describes a 
universal property of spacetime. Its knowledge allows, in principle, to deduce 
the constitutive relation of a material medium from the spacetime relation and 
the (charge) structure of the material, see the appendix [Bl Finally, it is also 
possible, when interested in material media, to abandon the deductive point of 
view and adopt a phenomenological description of the constitutive relation be- 
tween the macroscopic field excitation related to the macroscopic (free) charges 
and the field strength without trying to deduce it from the underlying structure 
of the system. From this point of view, the only physical property which differ- 
entiates what we call vacuum and what we refer to as a material medium, is that 
the former does not (by definition) have any charge substructure. It is a truly 
physical question whether the charges described by the charge density J are 
really the total charges within a certain region, or whether they are macroscop- 
ically averaged charges. It is important to realize that in any case, as soon as 
the considered charge is conserved, a corresponding electromagnetic excitation 
H can be introduced by mean of the inhomogeneous Maxwell equations. As a 
result, H will be just as 'microscopic' as the corresponding charge current J is. 

In general, spacetime/constitutive relations H = H(F) can have any func- 
tional dependence, i.e, they can be non-linear and non-local. Non-local consti- 
tutive relations are common for material media because in general the macro- 
scopic electromagnetic properties at some point of the material are influenced 
by the field and charge configuration at other points of it. Additionally, the fi- 
nite propagation speed of the microscopical electromagnetic interaction between 
the different parts of the material leads to time-like non-locality. In many cases 
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(typically, for slowly varying fields), however, this non-local effects are negligible 
and the medium can be treated as if the field excitation H were depending only 
on the value of the field strength at the same point, i.e. with a local constitutive 
law. In what follows, we restrict ourselves only to this later case. 



2.6.1 Linear constitutive relations 

We now concentrate on the particular case in which the field strength and 
excitations are proportional. A great number of material media are known 
in which this is a good description for a variety of conditions, see for instance 
[35]. We also expect the spacetime relation (i.e. the vacuum relation) to be 
simple and, in particular, linear. Linearity is taken for granted in all traditional 
approaches to vacuum electrodynamics. 

Given a local frame i) a , with a — 0, 1, 2, 3, we can decompose the exterior 
forms H and F as 

H=±H a/3 ti° t A^, F=^F a ^ a A< (69) 



and write a general linear spacetime relation (constitutive law) as 

H *0 = \ K a / F is = \ e a0 -ys X l5ef> F ee , (70) 



where i a /3jS is the Levi-Civita symbol with egjjg := 1 and -^ a ^ s is called the 
constitutive tensor density, which is an untwisted tensor density of weight +1, 
carrying dimension [\] = [H] / [F] = q 2 /h or, equivalently, dimension of con- 
ductance or I /(resistance). 

From its definition, the constitutive tensor satisfies the following symmetry 
properties 

which means that it has 36 independent components. Due to these symme- 
try properties, the constitutive tensor can also be represented by 6 x 6 ma- 
trix, say x IJ i where each index I, J, . . . corresponds to a pair of antisym- 
metrized 4-dimensional indices, i.e., I — > [21*2], etc. In this notation, e.g., 
e IJ — » e w2jij2 _ One can enumerate the 6-dimensional indices according to 
J, J, ... = 01, 02, 03, 23, 31, 12. 

We want to identify some irreducible components of the constitutive tensor. 
Usually, irreducible pieces of a covariant object have different physical proper- 
ties. Remember that so far no metric structure has been introduced. However, 
we always have the Levi-Civita tensor density e a ^ 1& at our disposal, and thus 
we can decompose the constitutive tensor according to, see |6Ij . 

^ap^S = (l) x a/3 7 <5 + (2)^Q/37<5 + (3)^a/3 7 <^ (J,}) 

where the different irreducible pieces ^X, an d are determined by the 
symmetry properties 

(l)^ Q /3 7 5 = (l) x 75a/3 ; (1) X W?7<5] = Qj (73) 
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(3)^0/375 _ (3) x [af3-tS]^ (3) x a/3 7 5 _ {3)^y5a0^ (75) 

or, explicitly, 

(3) x a/3 7 5 ._ x [a/37i] 5 ( 76 ) 

(2)^/375 ._ 1 |" x 75a/3 _ x Q / 3 7'5] ; (77) 

(1)^/376 . = x a/3 7 <5 _ (2) x Q/3 7 5 _ (3)^0^6] _ (73) 

One can introduce an alternative, but equivalent, parametrization of the 15 
independent components of the asymmetric piece ' 2 'x m terms of a traceless 
second rank tensor S£ (thus, also with 15 independent components) [55] as 

(2)^0/97* = ga/Jeh^fl _ e 7ie[a | g e ffl j = 0. (79) 

Now, contracting (|T9"|) with the Levi-Civita symbol, one finds that 

Sf = \ e a ,se (2) X^, (80) 

which shows that the traceless tensor Sj 3 is uniquely determined by (|79|) . 
The third piece can always be written as 

Wjf&r' = a(x) e a ^ s , (81) 

where a(x) is uniquely determined to be 

a := y aPl sX aM& - (82) 

It is called an abelian axion field. 

Denote the contribution of each piece of the constitutive tensor to the exci- 
tation as Wjf, (2)# an d ( 3 )# ; suc h that 

W H al3 = ±Z alh s a) X lSee F tB> (83) 
and analogously for the other pieces. Then 

H=<VH+WH + <®H, (84) 

and 

^H = aF. (85) 

Using (O and j39]), we find 

dH = d^H + ^H + ^H) 

= d^H + ^H) +d{aAF) 

= d( (1) H + (2) H)+daf\F, (86) 
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showing that the axion piece contributes to the Maxwell equations only if da ^ 
0. 

Equation (|72[) represents the irreducible decomposition of the corresponding 
6x6 matrix with respect to the linear group into a symmetric traceless piece, 
an antisymmetric piece, and a trace piece. The Levi-Civita symbol serves as 
a kind of 'metric' in the 6-dimensional formulation (e y *" — > e IJ ) which is used 
to construct the traces of \ IJ . The constitutive tensor is therefore reduced as 
X = X+ X m 36 = 20+15 + 1 independent components, respectively. 

No further decomposition of the constitutive tensor is possible at this point, 
since no additional geometric objects are available. 

For later application, see section I3TTI we abbreviate the action of the consti- 
tutive tensor by defining the operator # : A 2 — > A 2 such that 

# 0:= \e afilS x' St0 W^, (87) 

for any 2-form Q with frame components £l a p. Then we can rewrite our space- 
time/constitutive relation (|7D|) as 

H = *F. (88) 

Each irreducible piece is expected to describe different aspects of the medium. 
Additional information about the different properties of each piece will be ob- 
tained from the study of the electromagnetic energy-momentum current and of 
wave propagation. Notice for example that, if the Maxwell equations 

cannot follow as Euler-Lagrange equations from a Lagrangian of the usual form 
V := H A F, since ( 2 'X 7^ drops out from V due to its symmetry proper- 
ties. However, our intention is to try to develop our electromagnetic theory 
as generally as possible. Therefore we want to also include systems for which 
no Lagrange density can be found. Typical examples of such kind of physical 
systems are those including dissipative effects. One can therefore expect the 
'extra' irreducible piece to be related to some kind of intrinsic dissipative 
property of the medium. As we will see in section [2. 71 this is indeed the case. 
Furthermore, constitutive laws (for matter) with ' 2 'x ^ (non- vanishing "ske- 
won fields" ) have been discussed by Nieves and Pal [17J EE] • They yield T- and 
P-violating terms in the field equations. 

Non-abelian axions were postulated for the first time by Peccei and Quinn 
[55] , see also [82l [83] . Abelian axions coupling to electromagnetism were first 
considered by Ni [44] and correspond to a non- vanishing piece ^X- There have 
been intensive experimental searches for axions, see [151 [SSI US] and references 
therein. To date, no evidence of such a field has been found. Constraints on the 
axion mass and coupling to the electromagnetic field have been obtained both 
from astrophysical observations as from laboratory experiments, see [151 1561 175] 
for details. However, as we will see later, the axion-like term does not enter into 
important quantities as the energy-momentum current and the Fresnel equation. 
The axion remains a serious candidate for a particle search in experimental 
high energy physics and is a candidate for cold dark matter. The discussion 
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about the possible existence of such field for some material medium has been 
rather controversial, see [801 HH El] and references therein. In this context the 
vanishing on the axion piece is referred to as the 'Post constraint' (PC), after 
the work of Post [55] . Sihvola and collaborators have correctly recognized that 
the axion piece is allowed by the basic structure of electrodynamics. However, 
if the medium is homogeneous, a is constant and then the corresponding term 
drops out completely from the Maxwell field equations. Therefore, a possible 
axion field can only be detected if it is inhomogeneous (da ^ 0), or by its effects 
on the boundary separating two homogeneous media with different axion fields 
each. Some theoretical work on the reflection and transmission properties of this 
kind of medium can be found in [72] and references therein. No clear example 
of a material medium with a nontrivial axion-like term in its constitutive law 
has been reported in the literature. 

Within a \~9 formalism, i.e., when a (symmetric) constitutive tensor \ is 
assumed for the description of the electromagnetic properties of spacetime in 
addition to a metric <?, astrophysical constraints on 'nonmetric' theories, i.e., in 
which the constitutive tensor differs from the one on a riemannian space, were 
considered in [45] . The particular case of deviations from the constitutive tensor 
of Minkowski vacuum is studied in [201 HH 132] • 

2.6.2 Three dimensional decomposition of the constitutive tensor 

In some particular applications in which a 3 + 1 decomposition is used, as for 
instance in nonrelativistic and/or noncovariant formulations of electrodynam- 
ics, see section 12.6.31 for a particular example, it is convenient to express the 
constitutive tensor (36 components) in terms of four 3x3 matrices A, B, C, and 
V (each with 9 independent components). We defin^j] them as follows: 

c\ ■■= \ i acd x cdQb , v b a := 1 x 0acd ic db - (90) 

Here a, b, c, . . . = 1, 2, 3. The inverse relations are 

X 0o0b = A ba , x° abc = T> d a e dbc , (91) 

abOc abdnc abed £ a ^ e K c f c d (QO\ 

X — e ^ d> X — t °/e e • l yz J 

Additionally, one can (3 + l)-decompose the coordinate components of exci- 
tation and field strength as 

V a :=(H 23 ,H 3U H 12 ), H a :=(H 01 ,H 02 ,H 03 ), (93) 

B a := (F 23 , F 31 , F 12 ) , E a :— (F w , F 20 , F 3Q ) . (94) 

8 by the very definition of an homogeneous medium. 

9 The conventions in these definitions are taken such that they are consistent with those of 

mi- 
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With these definitions, the spacetime/constitutive relation can be written as 

V a ._ _ A ba Eb + v a B b^ ^ . = _<* ^ + g^gb (Q5) 

From these relations, one sees that the components of A correspond to the com- 
ponents of the usual dielectric tensor e, whereas the components of B correspond 
to those of the inverse permeability tensor fi, see ([9]). Furthermore, using the 
6-dimensional notation, we have 

H I = (H a ,V a ) : F! = (-E a ,B a ), (96) 

and the spacetime relation is written as 

Hi = e IJX JK F K , (97) 

with 

13 _ ( Bab V a h \ rj _ r _ / 3 1 3 



* "I K A* J' e = e -=l Z £ )• ^ 



Then the irreducible decomposition of \ can be found to be 

B(.6) \ {V>a b + Pa 



C V J =( i, 2 ^ a A (ab{ a> ), (99) 



(2) 13 _ I B [ab] k{C b a -Vj>) . 

X -I l(p b a_ c a b) ^[afc] J- 11001 
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(V J = g(C c c + 2V)e". (101) 
We introduced the 3-dimensional traceless quantities 

n ■■= C\ - ~ C c c 6 a b , Tp b := 2? a b - i 2? c c <t (102) 

so that 

r a = 0, ?> Q a = 0. (103) 

In terms of S, see (j5U|). the second irreducible piece of the 3-dimensional 
matrices (l89 |) - (|90|) are found to be given by 



Therefore 



^C b a = -S a b + S b a S c c , ^V b a = S b a -5 a b S c c . (105) 



(2) u _ , WB ab (2) V a b \_( e abc S c +S a b -5 b S c c 

X "I (2) C a (2) A ab -\_ S a + S a S c _ e a4cgr ) l^bj 



b • Al J \ u b ~ u b u c 



and 



(2) v a = e *bc Sb Ec + ( Sb a _ d a Sc C- )B b^ 

{2) U a = e abc S b B c + (S a b ~S b a S c c )E b . (108) 
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2.6.3 The asymmetric constitutive tensor of Nieves and Pal 

Consider the particular case in which is the traceless part of the product of 
a covector uj and a vector v, i.e., 

= - i (w fc « fc ) (109) 



Then we have 



(2) A ba = v e abc Uc) (2)^ = _^ ^ ^ (UQ) 



(2) c b = _ 



U a V b + ~Sl (oJ V° - IV C V C 



(111) 



&V b a = uj b v a + ^ (cj v° - uj c v c ) . (112) 

As we mentioned, Nieves and Pal [13 05] discussed P- and T- violating gener- 
alizations of the Maxwell-Lorentz equations for isotropic material media. The 
constitutive tensor they study has a nonvanishing skew-symmetric piece ^'X- 
A 4-dimensional formulation of the constitutive tensor (Nieves and Pal used 
the 3+1 form of Maxwell's equations in cartesian coordinates) can be given by 
considering, in addition to the Minkowski metric r), a time-like vector v, which 
explicitly breaks Lorentz invariance. 

Consider, in addition to the usual 'metric' piece 

d) x «« = ,f^(r, ik rf l - rf k rf l ) = 2.[^-^ k rf\ (113) 
V Mo V Mo 

an antisymmetric piece ( 2 'x defined by (|79|) and (| 1 09|) for the particular case in 
which uj and v are, in cartesian coordinates uj a = (u>, 0, 0, 0) and v a = (v, 0, 0, 0). 
This choice implies in particular that the material will look spatially isotropic 
in this frame. The corresponding 3-dimensional constitutive matrices A ab , B ab , 
C a b , and T> a b , according to their definitions (|89[) and (|9H)) . are then found to be 

A -b 4 -e Q 5 a \ B ab = tfSab, (114) 
C\ = -^8t V b ^ = ^f8l (115) 

Here we have also used 7700 = c 2 , rj a b = —Sab, V°° = ° 2 an d r l ab = —S ah . 

With these constitutive matrices, the 3+1 decomposition of the constitutive 
law takes the form 

2 2 
V a ± e 8 ab E b + B\ H a = E a + ^ l 5 ah B b , (116) 



or, in the vector notation used by Nieves and Pal [471 EH] j 

D = e E + 7 B, H = 7 E + /i - 1 B. (117) 

Therefore, the constitutive relation of Nieves and Pal can be recovered provided 
7 = ujvc 2 /2. 

Since we have introduced an additional timclike vector, it is natural that the 
material will violate T invariance, since v defines a preferred 'time direction'. 
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2.6.4 Nonlinear constitutive relations 



We are mainly interested in linear electrodynamics. However, the formalism 
developed here can also be applied to nonlinear models, in the case that one 
concentrates on the properties of small perturbations of the electromagnetic con- 
figuration. This is, for instance, the case when one wishes to study propagation 
of electromagnetic waves in a nonlinear model. Examples of such models are the 
effective Heisenberg-Euler nonlinear theory 29 ,30,60 and Born-Infeld electro- 
dynamics [5]. In this cases, our linear formalism can be applied by defining an 
effective constitutive tensor. 

Consider an arbitrary local relation H — H(F). The corresponding Maxwell 
equation in absence of charges is written, in components, as 

dj {e^Hki) = 0. (118) 

Consider now the properties of a small perturbation AF of the electromagnetic 
field around some background configuration F. We write the total electromag- 
netic field strength F as F = F + AF. Then the field excitation can be written, 
to first order in the perturbation AF , as 



H ij (F) = H ij (F) > 1 9Hij 



2 dF kl 



AF kl . (119) 



F 



Inserting (|119[) into (|118[) and assuming that the background field F is a solution 
of pi8p . i.e. dj \e l i kl H kl {F)\ — 0, we obtain an equation for the perturbation: 

1 » u dH kl 



dj f -e ijkl 



, - / 0- (120) 

We can write this equation in the same form as in the linear case, i.e. as 

d 3 ( X % kl AF kl )=0, (121) 
with the 'effective constitutive tensor' 



ijkl \ ijmn 9H mn 

XcS ■ 2 dF kl 



(122) 



The tensor will, in general, depend on the local constitutive law and on 
the background field F. This result shows that most of the results obtained for 
linear constitutive/spacetime relations can also be applied to every local elec- 
tromagnetic theory provided one considers perturbations on some background 
solution. This is the case, for instance, of the propagation of waves in nonlinear 
media. 



2.7 Symmetries and Energy-momentum 
2.7.1 Symmetry of a linear medium 

Here we define the concept of a symmetry of an electromagnetic medium. This 
definition applies to the linear case, when a constitutive tensor \ is available. 
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To motivate our definition, consider the particular case in which the field 
configuration is such that the Lie-derivative of the electromagnetic field strength 
along some vector field £ vanishes, i.e. 



£ 5 F = 0. 



(123) 



The geometric interpretation of the condition (|123p is clear. It means that 
the electromagnetic field strength F is 'constant' along the direction £. More 
precisely, the integral J n F on some 2-dimensional surface Q2 (which is the 
natural physical quantity associated to the 2-form F) is invariant under the 
displacement of the integration region SI2 induced by £. This property follows 
from the very definition of the Lie derivative, see, for instance, the discussion 
in [59] and figure [4] for *B = F and n = £. 

Now, in linear electrodynamics the excitation is determined by the field 
strength F and the constitutive tensor x- It is clear that even if (the integral 
of) the field strength F is constant under the displacements defined by £, i.e., 
if (| 1 23|) holds, the excitation will not be constant unless the medium itself, i.e. 
the constitutive tensor, satisfies some condition. We will call this condition a 
symmetry condition for the medium. Clearly, the condition we are referring to 
is the vanishing of the Lie derivative of the constitutive tensor, since in this 
case (|123l) implies that H is also constant, i.e. C^H = 0. This lead us to the 
following definition: 

Definition: A linear electromagnetic medium is said to have a symmetry 
under the displacement induced by a vector field £ if the Lie derivative along £ 
of its constitutive tensor vanishes, i.e. if 



for some vector field £. 

Notice that this condition implies the independent vanishing of the Lie 
derivative of each irreducible piece of the constitutive tensor, see the general 
decomposition formula (|72|) . From this definition and the property about the 
noncommutativity of the Lie derivative, it follows that if £1 and £2 are two vec- 
tor fields describing symmetries of the medium, then the new vector [£1,62] is 
also a symmetry of the medium. 

This immediately raises questions analogous to those one encounters in the 
study of isometrics in GR, namely a) What is the maximum number n max of 
symmetries that a constitutive tensor allows? And then, assuming that n max 
is finite, b) What is the form of a 'maximally symmetric' constitutive tensor?, 
i.e., of a constitutive tensor allowing the maximum number of symmetries? 

It can be speculated that the maximally symmetric constitutive tensor can 
correspond to x — X{r]}> i- e -j to the (conformal) Minkowski vacuum, since we 
know that it is a highly symmetric case, see also the particular case below. 
This result would be interesting since it would provide a further way to 'de- 
rive' the vacuum spacetime relation by postulating that it has to be maximally 
symmetric. These questions will be investigated in the future. 




= 0, 



(124) 
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2.7.2 Riemannian case 



Here I show that in the special case for which \ = X{g}i i- e - f° r the 'vacuum' 
on a riemannian spaced, the equation (|124[) defining the symmetries of the 
medium is equivalent to the conformal Killing equation for the metric g, i.e. to 

£ e (|9| V V)=0, (125) 

or, explicitly 

9 mj d m e + 9 im dU j + rd m 9 ij = \g l] [2d m r + r9kid m g kl ] ■ (126) 



One proves this as follows. From (J223J) one finds 

c,xtl = 2(m 1/4 ^)a(M 1/4 ^ 

+2 (l^V) A (lff| l/ V fc ) - (127) 
Contracting (| 12T[) with \ \g\ gjij, one obtains 

(\9\- 1/4 9j l) ^Xg? =9 ik 9il£ 6 (l^V) +2£ C (l^V) , (128) 
so that the condition (|124p implies 

5^/^ (ifff'V) +2£ ? (|<?| 1/4 5 lfc ) =0. (129) 

Contracting this equation with gn~ one finds gjiC^ (ig^^ g^ = which, when 

substituted back into (|129p . results in the conformal Killing equation for the 
metric g, namely (|125p . This means that the conformal Killing equation is 
a necessary condition for £ to be a symmetry of the medium. On the other 
hand, from (|127p one sees that if £ is a conformal Killing vector then (|125p is 
automatically satisfied. Therefore, (|125[) is also a sufficient condition. 

In other words, in the riemannian case, all symmetries of the medium are 
conformal symmetries of the metric. 



2.7.3 Conservation of energy-momentum 

In physics conservation laws play a central role. Among other properties, they 
allow the definition of conserved quantities in terms of which the description 
of the system and its evolution becomes simpler. In particular, energy and 
momentum are quantities associated to any field. In GR and similar theories 
they are the source of the gravitational field. 

In field theory, the energy-momentum distribution and its flow are described 
by a covector valued 3-form which, in the electromagnetic case we will denote 

10 or a material medium in Minkowski space with an 'effective optical metric' g. 
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as S a . Usually, one also needs an additional vector field £ for constructing a 
frame-independent 3-form Q := E a £ a following the general scheme 'conserved 
quantity' ~ 'momentum' x 'vector field'. The vector field is usually related 
to some symmetry of the system. As an example of this general scheme recall 
the case of GR. There, a conserved quantity along the trajectory of a free test 
particle (i.e., a particle moving along a geodesic) can be defined provided the 
spacetime admits a symmetry. This symmetry is described by a Killing vector 
field £kj i-e. an isometry of the spacetime metric, so that L^9ij = 0. The 
corresponding conserved quantity is then given by the projection (contraction) 
of the momentum\^\ of the test particle of mass to, pi :— mgijdx^ /dr, along the 
direction of the Killing vector i.e. q :— pi 

Back to our original problem, the main property that Q has to fulfill in 
order to be interpreted as energy-momentum is that it has to be related to the 
corresponding force law through a derivative of the form £ Q / Q ~ dQ. If, for 
instance, £ happens to be a 4-velocity field, S, a f a can be interpreted as the rate 
of energy transfer to the particles via the Lorentz force. A second condition 
on Q is that it must be conserved under some conditions, i.e. dQ — 0. One 
expects again the condition for energy-momentum conservation to be related 
to the symmetries of the system, at least when a Lagrangian is available, as 
we know from the Noether theorem. In linear electrodynamics this latter can 
happen when the irreducible piece ^ 2 'X °f the constitutive tensor vanishes. 

Therefore, to find an adequate energy-momentum E a for the electromagnetic 
field, we try to express £ Q / a as a total derivative of some 3-form dQ, with 
Q ■= S a f. 

First, we use the expression for the Lorentz force law (|37|) . replace the current 
J from the inhomogeneous Maxwell equation (|35|) and 'partially integrate'. One 
finds 

T/a = C{e a \F)AJ 

= (£1*0 a J 

= (H\F)AdH 

= -dmF)AH}+HAd(t\F) 

= -d [(£J F) A H] + H A (£ e F) . (130) 

In the last step we have used the definition of the Lie derivative of F and the 
homogeneous Maxwell equation, dF = 0. 

The choice of Q that equation (|130| suggests, namely Q = — (£JF) A H, is 
not what we are looking for since it would lead to conservation of Q only if the 
field strength is constant in the direction £, i.e. if C^F = 0, so that the last 
term in (|130[) vanish. But this is clearly a too restrictive condition. 

To improve the situation we use now the fact that (|130|) can be rewritten as 

C a fa =d[Z-{i\F)AH] + H A (C 6 F) - dZ, (131) 
where Z is an arbitrary 3-form. 

11 Momentum is always a 1-form, i.e., a covector, see, for instance, figure 2 in |59| . 
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To be sure that one is covering all the available possibilities, one should 
consider the most general 3-form Z that can be constructed using the available 
objects. In our case, they are only F, H, and £ Under these conditions. 
Z can only be a linear combination of the terms (£JF) A H and (£,\H) A F. 
Therefore, we make the ansatz 

Z = ai ^\F)AH + a 2 {f\H)AF, (132) 

with arbitrary constants a\ and a 2 . This implies 

dZ = ai d(£\F) AH-a 1 (Z\F)AdH + a 2 d(€\H)AF-a 2 (t\H)AdF 

= ai d (£J F) A H - 01 (f J F) A J + a 2 d{£\H) A F 

= a x (£ f F) AH - ai eU + a 2 [C 6 H - (f\ dH)} A F 

= 01 (£fF) AH — ai CU + a 2 (L 6 H) A F - a 2 (£J J) A F 

= aiOQF) AH + a 2 {C i H) AF+(a 2 -ai)f / Q . (133) 

Thus, we rewrite (|130l) as 

r/a = d[(a 1 -l)^\F)AH + a 2 {i\H)AF] + (l-a 1 )HA(C f F) 

-as A F - (a 2 - ai) f / Q . (134) 

This expression, which is an identity as soon as F and satisfy the Maxwell 
equations and f a is given by (|37p . leads us to define 

Q ■= z [(oi — 1) (CJF) A H + a 2 (CJ H ) A F] , (135) 

1 — ttl + CL2 

X := [(1 - ai )HA (£ t F) - a 2 {L t H) A F] , (136) 

1 — ai + a 2 

so that 

£ a f* = dQ + X. (137) 

The identity (|137p summarizes all the possibilities to write £ Q / Q as a total 
derivative of the available fields plus some rest. The problem now is to find 
suitable values for a% and a 2 so that Q and X satisfy our requirements. In 
particular, we need X to vanish under some reasonable circumstances. 

Up to here our results in this section are valid for any electromagnetic 
medium, since the spacetime/constitutive relation has not been used. To make 
up our minds about a reasonable choice of the constants it is useful to study 
the particular case of linear electrodynamics and compute X explicitly in terms 
of the field strength and the irreducible pieces of the constitutive tensor. 

We use C^F = \ C^Fij dx % A dx^ and similarly for H, C^H = | C^H^ dx % A 
dx 3 . Furthermore, we introduce the abbreviation a := 1/[4(1 — a\ + a 2 )] and 

12 In principle, one could also use the potential A as a further object. However, it seems 
that no 3-form linear in £, as required by the left hand side of 113H . can be constructed using 
this 1-form together with F and H. 
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compute 

X = a[{l-a 1 )H ij {C^F kl )-a 2 {C^H i j)F kl \dx i A oV A dx k A dx l 
= a [(1 - a x ) Hit {^ F ki) ~ ^ (A^tf) e ^ 
= 2a [(1 - ai) W fei (C^F kl ) - a 2 (C(H kl ) F kl ] e 
= a [(1 - oi) X Hij F ij (CtFu) - a 2 (C a kl ^ F tJ ) F kl ] e 
= a [(1 - oi) X^Fij (JC $ F U ) - a 2 fax"**) *ijFu ~ *iX m (C^) F kl ] 
= a [[(1 - oi) X kUj ~ a 2X ljU ] Fij (L e F kl ) - a 2 (L £X kHj ) F l0 F kl ] e 
"(1 - ^ - a 2 ) ^x Mij + {3) X kHj ) F l3 {C^F kl ) 
+ (l-a 1 +a 2 )^x kUi F ij (£ i F kl ) 



-a 2 £ 5 ( (1) x kUj + i3) X klij )F ij F M 



(138) 



Therefore, we find 



X = — Fir. 



(1 - ax - ag) /(x 
(1 - ax + a 2 ) 

a-2 



(1 - ax + a 2 ) 



( (1) X feW + (3) X fcHi ) -Ffei e. (139) 



From this result, one recognizes first that no choice of a\ and a 2 makes X to 
vanish in general. Thus, no conservation law is possible, unless additional con- 
ditions are fulfilled. This is, however, what one expects. One also sees that the 
irreducible pieces enter differently in (| 1 39[) . In particular, notice that the fac- 
tor (1 — ax + a 2 ) in front of the irreducible piece ^ X cancels completely. This 
means that the term ( 2 )^ fc ^ Fij (C^F k i) will always be present in the decompo- 
sition of X, no matter which values of a\ and a 2 one chooses. In other words, 
the irreducible piece ^ 2 ' X of the constitutive tensor will always prevent Q from 
being a conserved quantity (in charge- free regions), unless of course a much 
more restrictive condition, like C^F = 0, is satisfied by the field configuration. 
This result is consistent with our interpretation of ^x a s related to intrinsic 
dissipative properties of the medium. 
The 'best' one can do is to choose 



1 — ax — a 2 = 0, 



(140) 



so that the terms proportional to the symmetric irreducible pieces and ^x 
in the first term on the right hand side of (|139[) vanish. In this case, after 
substituting a 2 = 1 — a\ into (|135[) and (|136p . one finds that Q and X do not 
depend on ax, so that one obtains a unique result, namely 



i£\F) A H] , 



(141) 



X=-[HA (£ € F) - (C^H) A F] . 



(142) 
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Furthermore, for linear electromagnetic media we obtain, from (|139[) . 



X 



.(2 



X ijkl (C^F kl )+£ c («x ijW + <V W ) 



(143) 



That the choice (|141l) and (|142|) is reasonable can be seen as follows. Con- 
sider the case in which ^x = so that a Lagrangian is available for the de- 
scription of the system. In a charge-free region, (|143|) shows that Q is conserved 
provided £ is a symmetry of the medium, i.e., CfX — 0- Furthermore, our 
choice agrees with the results valid in a Riemannian space. 

Notice also that Q does not depend on the axion-like piece ^x entering the 
linear constitutive law. This can directly be seen using (|141|) and {85]), since 



= fcj (aF))AF-(Z\F)/\(aF) 

= a (£J F) A F — a (£J F) A F 

= 0. (144) 



However, the axion-like piece does appear in (|143p and therefore can contribute 
to the non-conservation of Q provided £^ a ^ 0. 



2.7.4 Axiom 4: Energy-momentum tensor 

The results in section 12.7.31 show how to construct a 3-form which for a linear 
electromagnetic medium leads to conservation laws, under reasonable condi- 
tions. We take these results and generalize them to our fourth axiom, by means 
of which the kinematic energy-momentum 3-form k E a of any electromagnetic 
medium is defined by 

k E Q := ~ [F A (e a \ H) - H A (e a \ F)\ . (145) 

Then the 3-form Q of section [2.7.31 is given by Q = k S a £". 

The adjective 'kinematic' is used here to emphasize that our definition of k S Q; 
was not based on dynamical properties of the system, but rather on kinematic 
arguments. 

One can also now define the components ^T a " of the kinetic energy-momentum 
tensor density by 

00 A k E a =: k Tj e, (146) 

where e :— $° A i) 1 A d 2 A d 3 . In particular, in a coordinate basis, with Ti.^ := 
e ijkl H kl /2, 

k V = ±6> i F kl H u -F lk W k . (147) 
For a general electromagnetic medium, we summarize our results as 

£ Q /a = dQ + X, (148) 

with 

Q:= k £ Q C*, (149) 



36 



X = 1 [H A - foff) A F] . (150) 

For a linear electromagnetic medium, we already showed that the axion-like 
piece ( 3 'x does not contributes to the kinematic energy-momentum tensor, i.e., 

k E Q ( X ) = k E Q (« X + + < 3 >x) - k S« ( (1) X + (2) x) • (151) 

Additionally, we have shown, see (|139[) and (| 148[) . that the corresponding 
conservation law takes the form 



dQ = f*C + I Fii 



A ( W X m + (i V' W ) F kl + 2^x ijM (£{Fh)] c, (152) 



for any vector field £. In this equation the contributions of the three irreducible 
pieces ^'Xi ^ 2 ^x, and ^ 3 ^x are isolated from each other. 

From these results, a consistent interpretation of (I152j) can be given as a 
'balance equation' for the energy-momentum content of the electromagnetic 
field. 

Recall first that even if £ is a symmetry of the medium, in which case the 
second and third terms on the right hand side of (| 1 52|) vanish, there is still a 
contribution proportional to the irreducible piece ^x °f l ne constitutive tensor. 
This is interpreted as an intrinsic dissipative property of the medium generated 
by the irreducible piece ( 2 >x- 

In the case that £ is a 4-velocity field the result above can be interpreted as 
an 'energy balance' equation for the change of the total electromagnetic energy. 
The integration of (| 1 52[) over a 4-dimensional region of the form 0,4 = £ x [ao, c] 
will produce at the left hand side the 'change of the total energy of the elec- 
tromagnetic field' between the 'times' ao and a, i.e., J E Q — J E Q. The first 
term on the right hand side, will lead to the energy transfered from the 

electromagnetic field to the test particles via the Lorentz force, i.e., J Q / a £ a - 
The term proportional to ^x can be then interpreted as the rate at which 
the energy of the electromagnetic field is dissipated, \ f u Fij^ \ l i kl (£{-Ffc;)e- 
The other terms will be proportional to the 'time' derivative of the constitutive 
tensor, | f n F^C^ (i x )^ kl + ( 3 )^y' fci ) F M e. When the medium is 'time' inde- 
pendent this term vanishes. If the time derivative does not vanishes, it means 
that the properties of the medium are changing in 'time'. Therefore, one could 
try to interpret this term as the energy per unit time needed to change the 
material properties 



3 Wave propagation 

We turn our attention to the wave propagation properties in our general frame- 
work of linear pre- metric electrodynamics. 

13 We assume that the fields vanish at 9E CT ('spatial infinity'). 

14 Maybe with a change of some suitably defined 'internal energy' of the medium. 
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Consider a region in spacctime without charges, i.e., J = 0. Maxwell's 
equations take then the form 

dF = dH = 0, (153) 

completed by the linear spacetime/constitutive relations (|70|) . These equations 
will allow for solutions propagating in spacctime, the behavior of which is de- 
termined by the electromagnetic properties of the spacetime/medium. 

Since the constitutive tensor can have a very complicated spacetime depen- 
dence, many specific features of the propagation of waves over finite distances 
cannot be studied in general terms here. However, an important local property 
of the propagation of waves, namely the dispersion relation that the covectors 
tangent to a wave front must satisfy, can be derived in general. This is done by 
deriving the so-called Fresnel equation for the wave covectors. 

3.1 Propagation of singularities 

In the theory of partial differential equations, the propagation of waves is de- 
scribed by Hadamard discontinuities of solutions across a characteristic (wave 
front) hypersurface S [T7]. One can locally define S by the equation Q(x l ) = 
const, for some function $(x). The Hadamard discontinuity of any function 
T{x) across the hypersurface S is defined as the difference {F] s (x) := F(x + ) — 
Fix-), where x± := lim (x ± e) are points on the opposite sides of S 3 x. An 

electromagnetic wave is described as a solution of the vacuum Maxwell equa- 
tions () 1 53(1 for which the derivatives of H and F have discontinuities across the 
wave front hypersurface S. 

Thus, in terms of field components, we have on the characteristic hypersur- 
face S, 

[Fij]s = 0, [diF jk ] s = qi f jk , (154) 
[Hij]s = 0, [diHjk\s = Qi hjk, (155) 

where and hij are the components of the 2-forms / and h describing the 
corresponding jumps of the derivatives of field strength and excitation across S, 
respectively, and the covector tangent to the wave front is given by 

q:=d<f> = qidx\ (156) 

Notice that (1154)) and (|155[) , taken together, are covariant conditions. In partic- 
ular, although (|154b ) and Q155b ) may not look covariant at first sight, they are, 
provided Q154h ) and Q155h ) are also satisfied. This can be seen by considering 
the transformation law of diFjk under coordinate transformations. Consider the 
coordinate transformation x l — > x l . Then we have 

d ( dxi dx k \ dx l dxi dx k 

di J ' k ' = cV 7 [dxl 7 dx^J jk + dx^dxTdx* diFjk - (157) 
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0(x)=const 





:[//]=[ x >=0 
*0, [dF] # 



Figure 7: Characteristic surface for propagation of electromagnetic distur- 
bances. The 1-form q := d<f> corresponds geometrically to two small parallel 
planes tangent to S at each point x. 



Assuming that the coordinate transformation is smooth across S, i.e. [-^7 



0. 



[d v Fj> k >] s 



s = 0, we find that the jump of diFjk across S is given by 



d ( dxi dx k 



dx 1 ' \dxi' dx k ' ) ^ jfe ^ s dx 1 ' dxi' dx k 



dx l dxi dx k 



(158) 



which reduces to the transformation law of a tensor field, provided (|154a ) is 
satisfied, i.e. 

dx 1 dxi dx k 



[di'Fj'k'] s = 



- [diF jh 



(159) 



dx 1 ' dxi' dx k ' 

We use (|154|) . (|155j) and the Maxwell equations (|153[) to find the compatibility 
conditions that the field discontinuities most satisfy, and find 



q/\h = 0, qAf = 0. 



(160) 



This ensures that the Maxwell equations (jl53[) are satisfied. The Hadamard 
method used here is equivalent to the usual geometric optics limit made by 



expanding a solution for the electromagnetic potential of the form A = ae 1 *. 
The covector q := d$ corresponds thus to the wave covector. 

Assuming now that the constitutive tensor is regular on S (in a geometric 
optics limit this means that we assume the scale of variations of the constitutive 
tensor to be much bigger that the scale of variation of the wave field) we obtain 
from (|70p the corresponding relation between h and /, namely 



h = *f. 



(161) 
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Substituting (|161[) into (|160a ) and using (|160b ). we see that the equations (|160[) 
reduce to 

<ZA # / = 0, 9 A/ = 0. (162) 
Now, the general solution of (|162b ^ is of the form 

/ = «Ao, (163) 

where the 1-form a is defined up to 'gauge transformations' 

a —>■ a + tpq, (164) 

for any scalar function tp. This is an algebraic consequence of the gauge free- 
dom in the definition of the electromagnetic potential. Then the compatibility 
condition (|162a ) reduces to 

qA*{qAa) = 0. (165) 

Furthermore, in order to allow for solutions, the fields / and h must satisfy 
some 'intcgrability' conditions. Since the form of the general solutions of (| 1 60[) 
is / = q A a and h = q A 6, we find 

/A/ = 0, hAh = Q, fAh = 0, (166) 

or, in components, 

e ijkl fijfkl=0, ev kl hijh k i = 0, e ijkl f ij h kl = 0. (167) 



3.2 Presnel equation 

As has been noticed before, see [5T], not all the equations in (|165|) are indepen- 
dent. This fact makes the derivation of the Fresnel equation more involved. 

Let us isolate the trivial parts of (|165p . We use a covector basis $ a in order 
to write the covector a in terms of its frame components a a as a = a a 'd a and 
choose the covector i?° as the covector q, i.e. # = q. This can, of course, 
always be done at each point of spacetime, since we assume q ^ 0. With this 
choice, a gauge transformation (|164[) can then be completely accounted for by a 
transformation — > ag + ip. This means that the other frame components cia, 
A, B, . . . = 1,2, 3, are gauge invariant. Since the term coming from the basis 
covector $° identically vanishes, i.e. gA*(gA i?°) = 0, (|165|) can be written as 

qA*(q Ati B )a B = 0. (168) 

Furthermore, only 3 of the 4 components in the 3-form (|168|) are non-trivial. 
Multiplying (fT68|) by i? 6 , one finds 

i? 6 AqA*(qAi9 B )a B = 0. (169) 
The remaining 3 components are 

tf A AqA*(qA$ B )a B =0. (170) 
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This equation can be written as 



eW AB a B = 0, (171) 

with a volume element e that we choose to be e :— d° A i? 1 A d 2 A i? 3 (it doesn't 
matter very much which 4-form one uses, since the equation is homogeneous) 
and a 3 x 3 matrix W AB . Using ((ST]) and (fTTO)) we compute 

d A Aqf\*{qf\-d B ) = A l? 6 A (e^sX 1 ** 6 &f &* A 
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= x mB i 

= - X 6AdB e. (172) 
Thus, since the negative sign is unimportant, we can define 

W AB. = X 0A0B (173) 

The corresponding Fresnel equation is then obtained from the vanishing of the 
determinant of W, as the necessary and sufficient condition for existence of non- 
vanishing solutions as- Now, since W is a 3 x 3 matrix, its determinant is given 

by 



W := det(W) 

\Ubc £def W AD W BE W Cf 



3! 

= ^ABctDEFX^X^X 6 ^- (174) 

We apply now the following procedure to rewrite W as a fully 4- dimensional co- 
variant expression. We 'complete' the 3-dimcnsional Levi-Civita symbols cabc 
to obtain the 4-dimcnsional one, by first using cabc = ^oabc an( l then taking 
one of the O-components of the constitutive tensors as fourth summation index. 
This leads us to consider the following identity: 

~ . pAQU. fiBOE. OCOF _ 1 £ n(30D ft~/OE tiSOF fT7K\ 

e OABCX X X =2 e "/37<5X X X • (175) 

This identity holds because on the right hand side of (|175|) , due to the properties 
of the Levi-Civita symbol, one of the indices a, (3, 7 or S must be zero, but on the 
other hand only a and will contribute, due to the (anti)symmetry properties of 
the constitutive tensor, see (|71|) . which would otherwise make one of the two last 
X-factors vanish. Finally, the two remaining contributions are equal, canceling 
the factor 1/2 and proving the identity. This allows us to rewrite (|174|) as 

vv - Y\2 ^ X X X ■ (17t>) 
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We repeat the same procedure to complete the remaining 3-dimensional Levi- 
Civita symbol. Now, we can use the following identity: 

WX"^X M V 56F = \ h P ,rX af3&P X^X 65XT - (177) 
Using ((177)) in p!76|) . we finally obtain 

W = -e^ h^x^X^X 6 ^- (178) 

Since d° = q = qi dx 1 , the above result can be written in coordinate components 

as 

W = ^ e mnM e rfltu x m "" X Jpsfc X' 9tn 9«&<fe«, (179) 

with := det(ef). We define the fourth order tensor density of weight +1, the 
'Fresnel tensor' Q % i kl , as 

. = i. g mnw e rstuX mnr(i X jlps]k X l)qtu , Q t]H = G im) , (180) 

which has 35 independent components. Then the Fresnel equation can finally 
be written as 

Q iM <nmm = o- (i8i) 

At each point, the Fresnel equation (|18ip defines in the space of wave covec- 
tors the wave ( co-)vector surface, see for instance [38]. I would like to emphasize 
the generality of the above result. The Fresnel equation (|18ip is valid for any 
linear electromagnetic medium. This means for media that in general are in- 
homogeneous, anisotropic, and dissipative. As we have already discussed, this 
result can even be applied to study the properties of propagation of pertur- 
bation of any local medium, by using an effective constitutive tensor, see |53j . 
Furthermore, the above result is, as the whole formalism, generally covariant. 

The general Fresnel equation (|18ip is in general a quartic equation in qi 
despite the fact that it was derived from a determinant of a 3 x 3 matrix quadratic 
in the wave covectors. This is because the remaining quartic term is multiplied 
by the trivial factor 2 , see (|179|) . This corrects Denisov & Denisov [14] who 
claim that a particular case of the general linear constitutive law may yield a 
sixth order Fresnel equation. 

3.2.1 Coordinate 3 + 1 decomposition of the Fresnel equation 

We rewrite now our general result (|179p by performing a 3 + 1 coordinate de- 
composition, and obtain 

W = 9 2 (q A M + q 3 q a M a + q 2 q a q b M ab + qoq a q b qc M abc + ladhdcdd M abcd ) , 

(182) 

with 

M:=g 0000 , M a :=4g 0Q0a , M ab := 6£ 00ab , (183) 
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j^j-abc ^QOabc j^j-abcd Qabcd 

or explicitly in terms of the 3x3 matrices defined in ([55)) and (|5D|): 



M = 
M a = 

M ab = 



(184) 

(185) 
(186) 



det„4, 

-hod (A ba A ce C d e + A ab A ec V d ) , 
i A^ [{C d d f + (V c c ) 2 - (C c d + V d c ){C d c + V c d )] 
+(C d c + V c d )(A c{a C b) d + v/ a A b)c ) - C d d A c{a C b) c 

_ Vc (a A b)c Vd d _ A dc c (a cV( b) + (A{ab) A dc _ A d(a A b)c^ g^g?) 



M a 



be 



de(c\ 



B df (A ab) Vj - V e a A b)f ) + B fd {A ab) C{ 



-A fla C b l) + C a f 2? e b) Vj + V f a C b \ C f d 



M 



abed 



e ef(c e \gh\d 



f "d ~r ^ e 1 - d 

- A ab "> K — C a T) b ) 



(188) 
(189) 



Each 3-dimensional tensor M is totally symmetric, i.e. M ab = M^ ab \ M abc = 
M (abc)^ M abcd = M (abcd)_ Thege results have been verine d by using the Maple 

computer algebra system, together with its tensor package GrTensor B 
3.2.2 Properties of the Fresnel tensor density 

First, one notices that the Fresnel equation is independent of the axion-like piece 
fi'X °f the constitutive tensor, 

0«tt( x )=g«*(d> x +( 2 ) x ), (190) 
which, in particular, implies 

g ijkl (^ X ) =0. (191) 
Furthermore, due to the antisymmetry property of ^Xt one verifies that also 

g ijkl (W X ) =0. (192) 

Actually, properties (|191|) and (|192|) generalize to 

g ijkl ( (2) X+ (3) x) = 0, (193) 

which can be verified, for instance, by using computer algebra. Notice that 
this identity is not trivial since g depends cubicly on the constitutive tensor 
X- The identity (|193|) shows that the symmetric piece ^'x is indispensable for 
obtaining well-behaved wave properties: If = 0, the Fresnel equation is 
trivially satisfied and thus no light cone structure could be induced. 



3 See http://grtcnsor.org 
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Furthermore, since 

g ijkl ( w x + (2) x) Q m ( {1) x), (194) 

the 'skewon' field does influences the Fresnel equation, and therefore, eventually, 
the light cone structure. An example of this general result can be found again in 
the asymmetric constitutive tensor studied by Nieves and Pal, see section [2.6.31 
and references [UJ BE] • 

Actually, after some algebra, one finds 

gW { W x + (2) x) = gW { {l) x) + | imnpg ^ (D x mnr(i (2)^.1* (2)^ 
+ \ e mnpq e rstu (2) % -nr(i (l) x i|p.|fc m x l)ltu^ (195) 

or, in a (more of less) obvious notation (see the definition (|180l) ). 

^ feZ (x,x,x) = g^ l ( {1) x, (1) x, {1) x) + 2G^ l ( {1) x, i2) x, {2) x) 

+^( (2) X, (1) X, (2) X). (196) 

The other terms vanish due to the symmetry properties of each irreducible piece. 

Take now (|f 95j) and substitute the parametrization of ( 2 'x m terms of S/ , 
see (f79|) . After some lengthy but straightforward algebra, one finds that the two 
last contributions to the right hand side of (|195[) are actually equal, namely 

S«*((i> X) eo X}( 2) x) = o««(< a >x, (1) x, (a) x) = iWx^i^sX . (197) 

Therefore, the final result reads 

5 iifcl ( (1) X + (2) x) = g ijkl { {l) x) + (1) x ra(i|n|i S m fe 5„ () , (198) 

a very simple expression, indeed. This equation summarizes how the skewon 
piece 'perturbates' the Fresnel equation that one would obtain only from the 
principal piece The skewon piece adds a second term which is linear in 

V-'X an d quadratic in ^x- 



2) 

For the particular case of Nieves and Pal, see ([79]) and (| 1 09[> . one finds 



G ljkl q t qjq k qi=-J— 
V Mo 



— (q ■ qf - (q ■ q){v ■ v)(v ■ q) 2 + (v ■ q) A 

Mo 



(199) 



Compare this result with equation (5.7) in [15]. Use qi = (uj,—k), v l 
(v, 0, 0, 0), gij = r/ij = (c 2 , —1, —1,-1), c 2 = l/eoMo and obtain 



[uj 2 - k 2 c 2 Y + ^k 2 uj~r 



(200) 



This results agreed 16 ! with that in [ISJ (in their e = /j = 1 case) for £ = cv 2 /sq. 



16 Nieves and Pal use a different system of units, compare, for instance, their vacuum equa- 
tions (2.1) and (2.2) with our Maxwell equations (Tf} and (|2}. Their equations can be obtained 
from our ones by substituting V = X>n, Ti = Wn/ c i E = E-^/eq, B = B-^/(ceo), p = p^, 
j = j'n, and t = tjf, where the subindex N refers to the fields in the notation of Nieves and 
Pal. 
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Notice that in general, for ( ^ 0, the Fresnel equation will have complex 
solutions. This is again a manifestation of the dispersive properties described 
by the antisymmetric piece °f the constitutive tensor. 

3.3 Light rays 

In a riemannian space, light rays are defined as the integral lines of the vector 
field defined by Vricm := g lJ Qjdi, which is then a null vector. It satisfies Vri em J q = 
0, so that the vector V T i em is orthogonal to the wave front S. This definition 
makes direct use of a metric, it is a nontrivial question whether it is possible to 
define light rays in our more general, pre-metric framework. 

In analogy to the riemannian case, one can try to define light rays as the 
integral lines of a vector field V which satisfies 

V\q = 0. (201) 

Equation (|201[) represents one condition for the 4 independent components of V. 
Therefore, this information is not enough to define light rays uniquely. We know, 
however, that in the riemannian case the light rays defined by V^em are also 
orthogonal to the polarization vector a (such that / = q A a), i.e. K-iomJa = 0. 
The same is true for b with h = q A b, i.e. Vri cm Jf> = 0. These known properties 
suggest to define the vector V for a given solution of the Maxwell equations as 
those satisfying 

V\q = 0, V\a = 0, V\b = 0. (202) 

The three conditions (|202j) determine then the vector V up to a scalar factor. 
This ambiguity in the definition of V is however irrelevant, since the integral 
lines of V are independent of that scalar factor 

One can construct immediately a solution of (|202[) . since they imply that 
the coordinate components V % of V have to be proportional to e v ' 1 qja^bi- We 
therefore define 

V 1 := e ijkl qj a k bi. (203) 

This quantity is a vector density of weight +1. Again, the fact that (|203p 
does not define a true vector field represents no problem for the definition of 
light rays as its integral lines. [Alternatively, one could try to 'normalize' V 
in order to construct a true vector, for example by dividing V by the density 
|a| 4 := aia^a^ai, with \a\ being a kind of 'norm' of a. Obviously, this 
procedure only makes sense if \a\ ^ 0]. One can write (|203[) in a coordinate- free 
way as 

V = ° (q A a A b) , (204) 

where ° is the (metric-free) operator mapping 3-form into vector densities of 
weight +1 according to ° (iT A ^ A ?9 7 ) = e Saf} ^ e s , so that V = V a e a . 

17 If one computes the integral line by dx 1 /dp = V 1 , then the freedom of choosing the scalar 
factor corresponds to the freedom of choosing the parameter p along the curve. In other words, 
an arbitrary scalar factor can always be absorbed by a reparametrization p' = p'(p). 
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Now, using the identities q A b = h = = &(q A a) we can write (|204p in 
terms of only q, a and the operator # as 

V = ° [a A * (a A <?)] , (205) 

or, in terms of the components of the constitutive tensor, 

V a = X af}lS apa iqs . (206) 

As a consequence of our definitions, V depends explicitly not only on the 
wave 1-form q but also on the polarization 1-form a. This property agrees with 
the result that, for a general constitutive tensor, different polarization states will 
propagate in different directions. Birefringence is again a particular example of 
the above mentioned fact. One can also verify that the axion-like piece of the 
constitutive tensor drops out from V, so it depends only on W% and ^x- 

To gain more physical insight on our definition of light rays, we evaluate the 
kinematic energy- momentum 3- form k Eo, for our particular solutions. From its 
definition (|f 45| and using F — ► / = q A a and H — > h = q A b, we find 

k S Q = i [q A a A (q a b- qb a ) - q A b A (q a a - qa a )} 

= —[q a q/\a/\b — q a q/\b/\a] 

= q a q A a A b. (207) 
Using (|204[) . we thus arrive at 

° ( k E a ) = q a V. (208) 

In terms of the components k T a @ of the kinetic energy- momentum density, see 
1111]), one finds 

V = q a vv. (209) 

This result shows that V determines the direction in spacetime in which the 
energy-momentum of the wave is transported. In other words, V represents the 
4-velocity of the energy transport. 

3.4 Presnel equation for V 

Now, since / = q A a and h = q A b, we see that the conditions (|202[) imply 

V\h = 0, VJ/ = 0. (210) 

Kiehn [34] uses (|210p to define V as an 'extremal vector field'. Furthermore, 
it can be shown that (|210|) are equivalent to (|202p . Take for instance (|2I0b ). 
then we have 

V\f = V\(qAa) = (V\ q)a-q (V\ a) = 0. (211) 

Suppose now that the scalar V\q does not vanish, then (|2 1 1 [) implies that a 
is proportional to q. But this would imply that / = q A a = 0. Therefore, 
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nontrivial solutions (/ ^ 0) require (|202a ) to hold. Substituting this into (|211[) 
shows that also (|202b ) must be satisfied. Similarly, starting from (|210a ). the 
same argument shows that (|202b ) is required. 

The equations (|210p for the vector V are similar to the Hadamard compati- 
bility conditions (|160p for the wave 1-form q. This suggests to apply a method 
similar to that used in section f3. 21 for deriving the Fresnel equation for q. 

First, we solve (1210b ) for /. A general solution can be written as 

f = V\c, (212) 

where c is an arbitrary 3-form. We write the 4 independent components of c as 
c Q( 3 7 =: iaf3jsc s - Using this in (|210a ) and (|88|) we find 

V\*(V\c)=0. (213) 

This equation is analogous to (| 1 65[) . As in the case of the wave 1-form q, not 
all the equations in (|2 1 3(1 are independent, since interior product of V with 
(|213p vanishes identically In analogy to the method used in section 13.21 it is 
convenient to use a vector frame basis e a and its dual i? a , with the special 
'gauge' ep = V. Then, after some algebra, one finds that (|213[) is equivalent to 

M AB c B = 0, (214) 

with the 3x3 matrix Mab defined as 

Mab := i 6 ACD X CDEF e EFM . (215) 

Again, nontrivial solutions (c B ^ 0) exist only if det (Mab) = 0, i.e. if and only 
if 

y€ ABC e DEF M ad MbeMcf = 0. (216) 

Using (|215p and (|216|) . and applying the same method as in section 13.21 to 
complete two Levi-Civita symbols, we finally find also a quartic equation that 
now the components of V must satisfy, namely 

M ijk i V l V ] V k V l = 0, (217) 

where 

Mijkl ■= ^ X mtef £efh(i£j\mnpX nPqr £qrs\k£l)tuvX UVhS > (218) 

is a totally symmetric tensor density of weight —1. 

Equation (|217[) defines in the space of ray vectors the ray surface, see |38) . 

One can verify after some algebra that in the riemannian case the tensor 
density M^ki reduces to 

M m := Mgf 1/2 g^gki)- (219) 
Additionally, one verifies from (|206|) that the vector density V reduces to 

V = A\g\ 1/2 [(g* k a k ) (g% qi ) - (g%) {g jk aj a k )} . (220) 

Since in a riemannian space we can additionally choose a 'Lorenz gauge' for 
a, namely g tJ aiqj = 0, the first term vanishes and the above result shows that 
V 1 ~ q l :— g^qj, as expected. 
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4 Light cone structure 



In this section, we want to study the particular case of our general theory in 
which a light cone structure is induced. In terms of the Fresnel tensor Q 1 this 
means that the spacetime/medium is such that 

sr* -\isP<r + M+0<t>), (22i) 

for some symmetric second order tensor density Q 13 of weight +1/2. Then the 
Fresnel equation, which in general is quartic in the wave-covcctor, reduces to a 
second order equation, namely to a light cone condition: 

Q t] q l q 3 = 0. (222) 

Obviously, we are interested in this subcase because in GR (and its special 
case of SR) the local properties of light propagation in vacuum are determined by 
the conformal structure of the underlying riemannian geometry. This implies 
that the wave covectors satisfy a condition ()222|) . with Q tJ ~ g ij , Here g 1 ^ 
denotes the spacetimc metric. 

We are specifically interested in the following questions: a) What is the most 
general constitutive tensor leading to a light cone structure? b) Is it possible 
to obtain a light cone structure as a consequence of some metric-independent 
conditions? In other words, can a light cone structure be induced/deduced 
without postulating the existence of a metric from the very beginning! c) If 
this is the case, how can a conformal metric be explicitly constructed from the 
underlying constitutive tensor? 



4.1 Looking for metric-independent conditions 

It is clear that in order to induce a light cone structure additional information 
or assumptions must the added to the general pre-metric framework, in the 
context of which the constitutive tensor is arbitrary. We look for these additional 
conditions. First, we consider the particular case of vacuum in a riemannian 
space in order to see if it is possible to find some properties which can be then 
used as a guide in the pre-metric framework. 

We know that the constitutive tensor corresponding to the vacuum in a 
riemannian space, i.e. 

^S} =2 ]ff ^S ilk 9 l]j , (223) 

is such that the Fresnel equation reduces to the light cone condition. Note 
that X{g} i s invariant under conformal transformations <?y — > e^^gij . As a 
consequence, only 9 of the 10 components of the metric can enter \- 
However, since the Fresnel equation is homogeneous, we see that 

V,;;;, ■>]■:,■■,—</',,'■. (224) 
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will also induce the same light cone structure, for any scalar function f(x). This 
field can be identified with a dilaton field. One uses the definition (| 1 80[) of the 
Fresnel tensor for this particular case and, after some algebra, we obtain 

<^ } =sgn( 3 )/ 3 /Ms ( V°. (225) 

Let us summarize the properties of the constitutive tensor (|224l) . Obviously 
^X{g,f} = 0, because X{g.f} is symmetric under the exchange of the first and 
last pair of indices. Furthermore, ^X{g,f} = 0, i.e., no axion-like term is 
present. Additionally, one verifies that 

1 2. ., , v klmn i v Mrs — f 2 f99fi1 
g £ ijkl X{gJ} tmnpq X{ g j} ~ J Z °i °j ■ 

In this expression the metric tensor does not explicitly appear, but the negative 
sign on the right hand side is due to the its lorentzian signature. 

These two properties, can be rewritten in terms of the corresponding linear 
operator k related to X{g./}i see (|70|) . The symmetry property ^ 2 'X{g,f} = is 
equivalent to 

A A (kB) — (kA) A B, (227) 
for any 2-form A and B, while the condition (|226[) can be written as 

k 2 = -fl. (228) 

From this example, which is realized in GR, we learn that the riemannian 
vacuum constitutive tensor satisfies some properties with can be written in 
a metric-independent way. As we will see in the next sections, an operator 
satisfying the two conditions above defines a duality operator. 

Hence the question arises whether fj22T[) and (|228[) can be taken as additional 
conditions in order to induce a light cone structure in the general pre-metric 
electrodynamic framework and whether these conditions are necessary and/or 
sufficient to induce a light cone structure. In fact, it turns out that there exists 
a relationship between dual operators and conformal metrics. 

4.2 Reciprocity and closure 

The two conditions (f227|) and (f228|) can be motivated already in the pre-metric 
framework. From the discussion of section 12. 7. 31 it is clear that the condition of 
symmetry ensures that the medium/spacetime will not possess intrinsic dissi- 
pative properties. This can then be taken as an additional metric-independent 
condition, see, however, the discussion in sections [2.6.31 and 13 . 2.21 about the con- 
stitutive tensor of Nieves and Pal, for which ^ 0. The condition of closure 
can be motivated as follows. An interesting feature of electrodynamics is its 
electric/magnetic reciprocity property. From the general expression of the ki- 
netic energy- momentum current (|145[) . one can see that the energy- momentum 
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content of the electromagnetic field remains the same if one exchanges excitation 
H and field strength F in an appropriate way. Consider the transformation 

H^(F, F^~H. (229) 

Here £ is a pseudoscalar carrying the same physical dimension as the constitutive 
tensor, i.e. [£] = [H/F] = q 2 /h. Direct inspection of (|145|) shows that the 
energy-momentum current remains invariant under this transformation, i.e. 

k E Q -► k S Q . (230) 

In terms of the 3-dimensional forms defined in section [2~4l the transformation 
H -> (F implies H -> -QE and V -> QB and F -> -± H implies E ^ ±H 
and -B — > — |Z> . This clearly shows that (|229[) exchange electric and magnetic 
quantities. 

So far, this invariance, which we call electric/magnetic reciprocity 0, is a 
property only of the energy-momentum current. If one requires, as an additional 
postulate, the spacetime/constitutive relation also to be electric/magnetic re- 
ciprocal, a condition on the constitutive tensor can be found. 

Take the spacetime/constitutive relation (1701) and perform the transforma- 
tion (I229p while leaving the constitutive tensor unchanged. One finds 

CF af3 = ^e a ^ sX ' tSee (-l)^H e e. (231) 

Use again (fTTTj) to replace the components of the excitation in terms of those of 
the field strength, and obtain 

- C 2 F a0 = ^ e aMS X lSte e^ X^ Xp F Xp . (232) 
This expression leads to a condition for the constitutive tensor, namely to 

~ C a *£$ = Jq W X^Ub^ X^ Xp - (233) 
Define now the dimensionless tensor 

£a/3 7 <5 ;= 1 X a/J 7 « ) ( 2 34) 

with 

C 2 := ~ {h 3 mn X mnm ) (e P9ra X raii ). (235) 
Then (|233p can be rewritten as 

\ e a ^6 X ^e^u X = -2*^ . (236) 



18 In order to distinguish it from the conventional metric-dependent duality transformation 
- - 7 — g ejjfcj g km g ln F mn , under which the vacuum (J = 0) Maxwell equations 



in a riemannian space are invariant 



50 



It is convenient to define 

J a p lS ■=\e aP r l A miS . (237) 
Then the closure relation (|236p reads 

\ J a ^ a J V a Xp = -2S[ X S P J (238) 

or, in the corresponding notation in terms of 6 x 6 matrices, 

J 2 = - 1 6 . (239) 

We call this condition on the constitutive tensor closure relation. The nega- 
tive sign in (|229b ) and therefore the one on the right hand side of (|238[) is a 
consequence of the negative relative sign of the two terms entering in the energy- 
momentum current, see (| 145[) . Below we will see that this negative sign will be 
responsible for the lorentzian signature of the induced metric. 

Mathematically, this means that the operator J represents an almost com- 
plex structure on the space of 2-forms. 

4.3 Dual operators and metrics 

A linear operator J : A 2 (X) — > A 2 (X) acting on 2-forms is said to be a dual 
operator defining a complex structure if it satisfies the property of symmetry, 

A A (JB) = (J A) A B, (240) 

for any 2-form A and B, and if it is additionally closed such that 

J 2 = -1. (241) 

Many people have studied this relationship, both as a useful tool in GR 
[6l [7] and as a method which could allow to consider the metric of spacetime 
as a secondary quantity constructed from some other fields, see for instance [8]. 
Clearly, the dimensionless part of the riemannian constitutive tensor / _1 X{g fc /} 
defines a dual operator which is just the Hodge dual of the metric g. 

In the context of electrodynamics, it seems that Peres [56] was the first to 
try to reconstruct the conformally invariant part of the spacetime metric from 
the excitation H and the field strength F. In his approach, however, the metric 
is assumed to exist such that the relation between H and F is just the same 
as the vacuum relation in a riemannian space. In our notation, this condition 
is equivalent to the postulate that the operator k is a dual operator and that 
it also equals the Hodge dual of some metric that is to be determined. In this 
sense, the approach of Peres can be considered as a sort of 'inverse problem'. 
Contrary to defining the operator k as the dual of the metric tensor, as is done 
in GR, Peres tried to determine a metric such that its Hodge dual operator 
coincides with k. To that goal, it is assumed that the operator k satisfies the 
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conditions of symmetry and closure, as a Hodge dual operator does. However, 
as we will discuss, Peres did not succeed in deriving the conformally metric since 
his result is not unique. 

Toupin [71] and Schonberg 62J also studied how a conformal metric structure 
is induced and in particular how the conformally invariant part of the metric 
can be deduced from the spacetime/constitutive relation under the assumption 
of symmetry and closure. They seem to be the first who were able to show that 
a conformal metric structure is actually implied as a consequence of symmetry 
and closure. Brans [5J [7] also recognized that, within general relativity, it is 
possible to recover the metric from its Hodge dual operator. These structures 
were subsequently discussed by numerous people, by 't Hooft [70], Harnett [IB] , 
and Obukhov & Tertychniy [50], amongst others, see also the references given 
there. 

From these studies, it seems clear by now that there is a one-to-one cor- 
respondence between linear operators (of the kind of k) satisfying symmetry 
and closure, and conformal metrics. If a metric is available, one can imme- 
diately construct the corresponding Hodge dual operator which, for a metric 
with lorentzian signature, will satisfy (|240p and (|241|) . On the other hand, the 
contrary is also valid. If an operator J satisfying (I240[) and (|24ip is available, 
then a conformal metric structure is induced, such that the operator J can be 
written as a Hodge dual operator. Furthermore, the metric components can be 
(re-) constructed by using the so-called Schonberg-Urbantke formula. Schonberg 
seems to be the first to derive this formula in the context of electrodynamics. 
This formula was also found by Urbantke [73l [74], but in a different context, 
namely in the framework of SU(2) Yang-Mills theories. We will discuss some 
derivations of this important formula below. 

Finally, by using the Schonberg-Urbantke formula in the context of linear 
pre-metric electrodynamics, Obukhov and Hehl [52] presented an explicit con- 
struction of the conformally invariant part of the metric tensor in terms of 
quantities parametrizing a constitutive tensor satisfying symmetry and closure. 

After discussing these developments, an alternative procedure to deduce the 
conformal metric and its light cone structure will be presented. It relies on 
the general results about the Fresnel equation governing the local properties 
of electromagnetic waves. This alternative approach is therefore much more 
suitable to understand the physics underlying the emergence of the conformal 
structure of spacetime, since it directly involves the properties of the propaga- 
tion of waves. Finally, this approach also allows to study how the dropping or 
modifying of the assumptions of symmetry and closure would affect the light 
cone structure. In particular, we will explore the consequences a possible asym- 
metric spacetime/constitutive relation (i.e., with skewon part) would have on 
the light cone, but maintaining the assumption of closure. 

4.3.1 Schonberg-Urbantke formula 

Here we discuss the successful derivation of of the conformal metric by Schonberg 
and the related work of Harnett. 
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In his work [55] Schonberg studied a different derivation of the metric from 
electromagnetism. Additionally, the analysis in [62] helps to better understand 
the group theoretical aspects of the emergence of a lorentzian metric structure 
and thus of the Lorentz group. 

As we already saw before, the electromagnetic field strength F, with its 
six independent components, can be described as a vector in a 6-dimensional 
real vector space. We call this space Sq. One can map a basis of the space 
of 2-form, dx l1 A dx l2 to a basis b 1 of Sq by means of the rule [1^2] — ► J = 
01, 02, 03, 23, 31, 12. Then one can write a vector A € S 6 as A = Ai b 1 . 

An important role is played by an inner product e on Sq induced by the 
4-dimcnsional Levi-Civita symbol e ufci . Given two vectors A and B of one 
can define their product as 

e(A,B):=e IJ A I B J . (242) 
In 4-dimcnsional notation, this is equivalent to 

e(A,B)e := AAB (243) 

or, in components, 

e(A,B) = ±e i i kl A ij B kl . (244) 

Thus, the Levi-Civita symbol acts as a metric in the 6-dimensional space Sq. 
Notice, however, that e(A, B) is a 4-dimensional density of weight +1. Since 
the eigenvalues of e IJ are +1 and —1, each with multiplicity 3 (see (|98b )). the 
signature of e IJ is (+1, +1, +1, — 1, — 1, — 1). This immediately shows that the 
6-dimensional space 5*6 naturally contains a SO(3, 3) group structure. Transfor- 
mation of S^-vectors under the action of the SO(3, 3) group leaves the product 
(|242p invariant. It is important to emphasize that this group structure is always 
present, independent of any metric, affine, or whatever additional structure on 
the 4-dimensional manifold. 

Using this product, one can express the assumption of symmetry of the 
constitutive tensor, i.e., ^\ — as 

e (A, JB) = e ( J A, B) (245) 

for all A, B e S6, because 

e (A, JB) = e IJ A! jf B K = X IK MB K (246) 
since x IK — J j K ■> ancl 

e (J A, B) - e IJ Jj K AKBj = X IK A K Bj = X KI MB K , (247) 

so that (I245|) requires symmetry of Xi i- e -i X IJ — X JI ■ 

Furthermore, the additional condition of closure of the operator J, i.e., 

J 2 = -1, (248) 
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is assumed to hold. 

In the 6-dimensional real space Sq, the introduction of the linear operator 
J satisfying (|245|) and (|248|) induces a almost complex structure. In order to be 
able to discuss this structure and in particular the eigenvectors of the operator 
J, consider the complex extension of 5*6, i.e. the vector space of the complex 

2- forms Sq(C). The negative sign on the right hand side of (|248|) implies that 
the eigenvalues of J are ±i, and since J is a real operator, each eigenvalue 
must have multiplicity 3. This means that the corresponding self- dual and 
antis elf- dual subspaces spanned by vectors with eigenvalues +i and —i are both 

3- dimcnsional. We denote these subspaces as S^(C) and S^(C) respectively. 
The vectors of (C) will be denoted as A^ 1 and satisfy 

JA ± =±iA ± . (249) 

Take any A ± £ S^(C) and write it as A* = A + iA' with A and A' real 
6-dimensional vectors (i.e. A, A' £ S$). Then from (|249[) one finds that 

J(A + iA')=±i(A + iA'), (250) 

with implies 

JA = tA', JA' = ±A. (251) 

Using (|25ip one can rewrite every element of S3 , in terms of only the real 
A £ S§ and the operator J, as 

A ± =AtUA. (252) 

The spaces Sq(C) and Sg(C) are orthogonal with respect to the product 
e. Take any A+ £ S£(C) and B~ £ S<r(C). Write them in the form ([2%2"j) . i.e. 
A + = A — i J A and B~ = B + iJB, and compute their e-product. Using the 
symmetry and closure properties (|245|) and (|248|) one finds 

e(A + ,B-) = e(A — iJ A, B + iJ B) 

= e(A,B)+ie(A,JB)-ie(J A,B) + e(J A, J B) 

= e(A, B) + ie(A, JB) - ie{A, J B) + e(A, J 2 B) 

= e{A,B)+e{A 1 -B) 

= 0. (253) 

Consider a real operator M acting on the real space 5*6. If A £ Sq then also 
MA £ Sq. However, the self-dual element (MA) + of will in general not be 
the result of the application of the linear operator M on A + , since 

(MA) + = {MA) - iJ(MA) ^ M(A - iJA) = M(A + ), (254) 

unless the operators M and J commute, i.e. [M, J] = 0. In other words, not 
every linear operator on Sq is a linear operator in S% . In particular, not every 
element of the group 0(3, 3) which leaves e invariant defines a linear operator 
in the self-dual space S3". 
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The corresponding subgroup of 0(3, 3) which commutes with J, i.e. the 
subgroup form by those operators N such that 

e (NA, NB) = e (A, B) , [N, J] = 0, (255) 

can be shown to be isomorphic to the Lorentz group. To prove this, consider 
the product e restricted to the self-dual subspace S% , which we denote as e + , 
and is defined by 

e+ {A+.B+) :=e{A+,B+) . (256) 

Consider operators N satisfying l|255p which are therefore also linear operators 
of Sg - . It is clear that they will also leave invariant the 3-dimensional product 
e + , since in this case 

e+ ({NA) + ,{NB)+) = e+ (NA+,NB+) =e(NA+,NB+) 

= e(A + ,B+)=e+ (A+ ,B+). (257) 

This shows that the subgroup of 0(3, 3) commuting with J forms the invariance 
group of the 3-dimensional metric induced by e + in Sg~. But since S$ is a com- 
plex linear space, the corresponding invariance group of e + is clearly SO(3, C), 
which is isomorphic to the Lorentz group SO(3, 1). In other words, an operator 
J satisfying symmetry ([IMS]) and closure induces a SO(3, C) w SO(3, 1) 

group structure, which manifest itself as the symmetry group of the natural 
metric structure on the self-dual space of J. Similarly, the same arguments can 
be repeated considering the anti-self-dual space . 

We consider now how to reconstruct the metric components of the corre- 
sponding induced lorentzian metric. Consider a basis for each subspace. We 
denote the basis of by fin? , a,b, . . . = 1,2, 3. Then the complex conjugate 

Si a) := (S^ )* is a basis of S£", see ([MS]) . Since S 6 (C) = S+ S3, the six 

Sg-vectors form a basis of Sq(C). Then the orthogonality prop- 

erty (1253P implies that the components of the metric e of Sq(C) in the basis 
js^, form a block-diagonal matrix, since e ^S 1 ^', S^^j = 0. 

One can use the exterior product of vectors in Sq(C), which we denote by 
A, and determine a volume element of S$ as 

W+ := A A = i l ahc S { f A Sf A sf . (258) 

Then W + is a 3- form on Sg, i.e. W + G A 3 (S , 6 ). Since S$ is 3-dimensional, 
W+ is, up to a scalar factor, independent of the choice of the basis S+\ Fur- 
thermore, any product of basis vectors Si" of an order higher than 3 vanishes 
identically (e.g., S { + } A S<_ 6) A A S { f = 0). A volume element for 
can be similarly defined by using the basis of S$ , which then satisfies 
W- = (W+)*. ' 
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A 3-form W G A 3 ^) has 20 independent components Wjjk, which corre- 
spond to a 4-dimensional tensor Wi 1 i 2 j 1 j 2 k 1 k 2 of order six, with the symmetries 

^iii2JiJ2fcife2 = ^Jijaiiia&ifei = ^^h\kijxjii\ 12 = Wjj.i2k1k2i1.J2 and, of COUrse, 
Wi l i 2 j l j. 2 k 1 k2 = 

— Wjaiiju'afcife j e t c - From this tensor, a symmetric second order 4-dimensional 
tensor density of weight +1 can be defined as 

Wir-=\e klmn W iklmnj , (259) 

which has thus 10 independent components. The remaining 10 components of 
W can be mapped to a second order contravariant tensor density of weight +2, 
defined by 

W ij :=\e lklm W klmnpq e n ^. (260) 
Following (|259[) we define the tensor density associated to W+ as 

W±:=\e hl ™W+ lmnj , (261) 

or, equivalently, 

H-,;: ^ 6 fei " m ^ abc S+W S+ (b) S+« . (262) 

Furthermore, we saw that a change of the basis of S% to a new one 

S'^ leads to a new volume WL which is necessarily proportional to W+, i.e. 
VF! = 03 W + for some (13. One can therefore define a new tensor density by 

H£:=(detW+r 1/4 W+ (263) 

Then W+ is a tensor density of weight —1/2 which is independent of the choice 
of the basis of S$ . Because of this important property, VV^ must describe 
an intrinsic property of the self-dual space and therefore of the operator 
J. It is then tempting to identify the symmetric tensor density (|263p with the 
conformally invariant part of the metric tensor |g| _1 ' 4 gy, which is also a tensor 
density of weight —1/2, i.e., 

W+=a| 3 r 1/4 ^, (264) 

with some (in general complex) factor a. 

That this identification is consistent with our expectations can be checked 
as follows. Consider the linear operator L acting on vectors of Se(C), defined 

by 

{LA) tj := l -L l3 kl A kh (265) 

with 

L i} U :=Wtm^ mnU - (266) 



56 



This operator is expected to be related to the original operator J, because the 
identification (|264[) would imply 

L t3 kl ^a 2 \g\- 1 / 2 g m g jn e mnkl (267) 

which is proportional to the Hodge dual operator of the metric g, when applied 
to 2-forms. 

Using (|263[) one re-writes the definition (|266p of the operator L as 

L t kl = (detW+)- 1/2 Wt m Wf n e mnkl . (268) 

Furthermore, it can be proved that L is a complex operator. From (|268[) and 
the definition (|262[) one can verify by direct but rather lengthy calculations that 
the elements of S$ and S$ are eigenvectors of L with eigenvalues +1 and — 1, 
respectively, i.e. 

LS + =S + , LS~ = -S~. (269) 

Since S% = (S£) , the operator L cannot be real. If it were, (|269a ) would imply 
LS~ = S~, contradicting (|269b h It is therefore clear that the real operator J 
is just given by J = iL, because 

(iL)S + = iS+, (iL)S~ = -iS~, (270) 

which coincides with (|249| . The property (|270|) is valid for any vector of S% 
and since Sq(C) = S^(BS^ , it implies that the action of iL and J on any vector 
of Sq(C) is exactly the same, hence they are the same operator. 
Thus, we have proved that 

J = iL, (271) 

so that, from (|268p . we have 

J % kl = (-detW+r 1/2 W+ m W+ n e mnkl . (272) 

With the identification (|264[) we write the operator J as 

V =V=99im9 jn e mnkl - (273) 

From this expression, one sees that the metric gij can be taken to be real, with 
lorentzian signature, since J is a real operator. In other words, the results and 
(|262[) allow us to reconstruct the metric components as 

klrnn ~ Q+( a ) Q+( b ) Q+( c ) (97 A\ 

with g^ real. The conformal factor remains of course undetermined. 
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4.3.2 Necessary and sufficient conditions for the constitutive tensor 

The results discussed above can be summarized as follows. If the operator J, 
defined in (|237[) by the dimensionless part of the constitutive tensor, is symmet- 
ric and defines an almost complex structure, i.e., (|240p and (j241[) are satisfied, 
then a lorentzian metric g is determined, up to a conformal factor, such that 

X l3kl = 2 v / ^/V li , (275) 

with g := det g^ < 0. In other words, the conditions (|240p and (|241|) are 

sufficient to be able to write X in the form (|275[) 

But it is clear, see the discussion in section 14.11 that (|275p is sufficient to 
ensure that the corresponding operator J satisfies (|240[) and (1241[) . We thus 
conclude: 

The necessary and sufficient conditions to be able to write the dimensionless 

o 

part X of the constitutive tensor as 

X t3kl =2^g l[k g l] \ (276) 
with a lorentzian metric g, are symmetry 

X ljkl =X kHj (277) 

and closure 

I iafijB X ^ 6 e e6tiv X ^ Xp = -28^5"} . (278) 



4.4 General solution of the Closure Relation 

In order to be able to find an explicit expression for the metric in terms of 
quantities describing the components of the constitutive tensor, we have to solve 
the closure relation (|233[) or, equivalently, (|238[) . We will solve this equations 
for the general case in which the constitutive tensor is asymmetric, so that we 
can later investigate the effect of relaxing the symmetry condition. 

Let us now make the closure relation explicit. We turn back to the constitu- 

o 

tive 6x6 matrix (|239p . We define dimensionless 3x3 matrices A '■= etc. 
In terms of these dimensionless matrices (we immediately drop the small circle 
for convenience), the closure relation reads, 

A ac B cb + C a c C c b = -61 (279) 

C a c A cb + A ac V c b = 0, (280) 

B ac C c b + V a c B cb = 0, (281) 

B ac A cb + V a c V c b = -S b . (282) 
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In 3 x 3 matrix notation, we then have 

AB + C 2 = -1 3 , (283) 

CA + AV = 0, (284) 

BC + VB = 0, (285) 

BA + V 2 = -1 3 . (286) 

Assume det2? ^ 0. Then we can find the general non-degenerate solution as 
follows: Define the matrix K a f, by 

K:=BC. : i.e. C = B^K, (287) 

and substitute it into (|285j) . then 

V = -KB- 1 . (288) 

Next, solve (|283p with respect to A: 

A= -B- 1 - B^KB^KB- 1 . (289) 

Multiply (|289p by C from the left and by T> from the right, respectively, and 
find with (j2"g?|) and ([2"ggj) . 

CA = -B^KB- 1 - B- 1 KB- 1 KB- 1 KB- 1 , (290) 

AV = +B- 1 KB- 1 +B- 1 KB- 1 KB- 1 KB-\ (291) 

Thus, (|284p is automatically satisfied. Accordingly, only (|286[) has still to be 
checked. We compute its first and second term of its left side, 

BA = -l a - KB^KB' 1 , (292) 

V 2 = KB~ l KB~ l , (293) 

and find that it is fulfilled, indeed. 

Summing up, we have derived the general solution of the closure relation 
(|239[) in terms of two arbitrary matrices B and K as 

A = -B- 1 -B^KB^KB- 1 , (294) 

C = B^K, (295) 

V = -KB- 1 , (296) 



or, in components, 



A ab = _j3 ab -B ac K cd B de K ef B fb , (297) 

C\ = B ac K cb , (299) 
V a b = -K ac B cb . (300) 
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Here B ab are the components of the inverse i.e. B ab B bc = 8 b - The solution 
thus has 2 x 9 = 18 independent components. Alternatively, one can write the 
solution in a more compact notation as 

A = -(la + C 2 )^- 1 , (301) 
V = -BCB' 1 , (302) 

which is parametrized by the arbitrary matrices B and C with altogether 18 
independent components. In components, this means 

A ab = -{5 a d +V c a V d c )B db , (303) 
B ab = B ab , (304) 
C b a = -B ac V d c B db , (305) 
V b a = V b a . (306) 

(307) 

4.4.1 Explicit derivation of the metric components from the consti- 
tutive tensor 

Here, the Schonberg-Urbantke formula will be used to find an explicit expression 
for the conformal metric in terms of quantities parametrizing a constitutive 
tensor satisfying the assumptions of symmetry and closure. From the symmetry 
condition one finds that the constitutive tensor has the form (l9"51). but with 



A ab = A ba , B ab = B bai V b a = C a b , (308) 

so that ( 2 >x — 0, see (|100P - These symmetry conditions restrict the general 
solution of the closure relation found in the last section. From (|299p and (|300[) 
we see that (|308b ) implies ~K bc B ca = B ac K cb . Multiplying this equation by 
Bad and using (1308b ) we obtain —Km = Kdb which tells us that the matrix K 
must be antisymmetric. Then (|297|) implies that A is automatically symmetric, 
so no further conditions on K arise from (|308b ). In summary, a symmetric 
constitutive tensor satisfying the closure relation is given by (|297() — (|300p with 

B ab = B ba , K ab = -K ba . (309) 

As we saw in section 14.3.11 to construct the metric using the Schonberg- 
Urbantke formula, we need a basis of the 3-dimensional space S$ of self-dual 
2-forms of the duality operator J. 

To that end, we decompose the basis b 1 of Sq into two 3-dimensional column 
vectors, according to 

b 1 = f ^ ^ , a, b,... = 1,2,3. (310) 
Then we can find their self-dual parts, 



b I + = ~(b I -Ub I ), (311) 
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and decompose them similarly as in (|310[) . into 



0% 
7 6 + 



(312) 



Assuming, as we did before to find a solution of the closure relation, that 
the matrix B is nonsingular, we can show that one can take 7 + as a basis of the 
self-dual space S$ , since the remaining self-dual 2-forms (3+ can be written as 
linear combinations of the former. This is expected since S3 is 3-dimensional. 
Explicitly, we have 

(3% = (i5 a b + B ac K cb )B bd 7 j- (313) 



This allows us to choose the 2-forms 7+ or, equivalently, the triplet 



(314) 



as basis of S% . 



The information of the constitutive tensor x is now encoded into the triplet 
of 2-forms S ( + } . 

The (Sl) zj are the components of the 2-forrn triplet = (<S+)^ dx i hdx i /2. 
If we substitute the self-dual 2-forms into (|274j) . we can display the metric 
explicitly in terms of the constitutive coefficients: 



9ij = 



detB 



- h 



ka 



B ab + {dct B)- l k a h 



(315) 



Here k a :— e abc K 0C /2, k a :— B a bk h and </> is an arbitrary factor. The determinant 



of this metric is found to be g = ■ 
in (|315[) has lorentzian signature. 



(det B) so that we verify that the metric 



4.4.2 Properties of the metric 

The inverse of (|3 1 5[) can be found to be 



1 / 1 - (det B)^k c k c 


-k b 


(j)detB\ -k a 


-{det B)B ab 



(316) 



With the help of (I315|) and (I316j) . we can compute the Hodge duality operator 
* attached to this metric. In terms of the components of the 2-form F, we have 



•p.. — 



9 ~ n^ m n^ n P 

tijkl y y £ mn • 



(317) 



This equation can be rewritten, in analogy to (|70|1 . by using the constitutive 
tensor X{g}i defined in (|223|) . so that 



* P. - — _2 



kill 



ijkl X{g} Fmn, 



(318) 
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In order to compare X{g} with the original constitutive tensor x, we compute 
the corresponding 3-dimensional constitutive matrices of X{g} according to (|89[) 
and O- 

Then straightforward calculations yield 

A t} = V^(g 00 g ab -g 0a g ob )=A ab , (319) 

(B {g} ) ab = ~^(g™g d f -g de g ef )e acd e efb = B ab , (320) 

(C {g} ) a b = 1^=9 {9° c 9 ad ~ 9 ac 9° d ) e bc d = B ad K db = C%. (321) 

Thus, X{g} — Xi i- e -i the metric extracted allows us to write the original duality 
operator J as Hodge duality operator, J = *, when applied on 2-forms. 



4.4.3 Alternative derivation 

Here we provide an alternative derivation of the conformal metric, this time 
based on a direct computation of the Fresncl equation for the symmetric solution 
of the closure relation given by (1297|) . (|300j) and ()309|) . Using this equation and 
our general expressions (I186j) - ()189j) for the Fresnel equation in terms of the 3- 
dimensional constitutive matrices, we find after a straightforward calculation 
that 

i A fc a fc ax2 



M = -detB^-detfiJ < (322) 

MCl = deb^^-cS)' ^ 

M ab = - — !— 4fc a fc b + 2B afc ( 1 - , (324) 

detB V detB/' v ; 

M abc = -AB b(a k c \ (325) 

M (abcd) = _^ et B )Q(ab B cd) ( 326 ) 



Substituting all this into the general Fresnel equation (jl82[) . we find 

W = 



% 2 (i-fefS]-2 9 o(gan 



detB ^ u V detB, 

-(detB) (q a q b B ab )] 2 . (327) 

Therefore we find that the Fresnel equation, W = 0, can be written as 

(ff y «ft<&) 2 = °. (328) 

where g*- 7 is the (inverse) 4-dimensional conformal metric which arises from the 
duality operator and the closure relation. Direct comparison with (|327p shows 
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that 



- K 1 "^)' < 329 > 

g 0a = -i>k a , (330) 

g ab = (detB) B ab . (331) 

Here tp in the undetermined conformal factor. 

Thus we indeed recover the null cone structure for the propagation of elec- 
tromagnetic waves from our general analysis. The quartic surface degenerates 
to the null cone for the conformal metric g. 



4.5 Relaxing the symmetry condition 

Having verified and explicitly constructed the metric tensor from the constitu- 
tive tensor in the case in which it is symmetric and satisfies the closure rela- 
tion, we want to study what consequences a nonvanishing skewon piece in the 
constitutive tensor could have on the emergence of a conformal structure. In 
particular, can the spacetime metric still be constructed? We have seen that 
the skewon piece does influence the light cone structure, and we also saw that 
an asymmetric constitutive tensor satisfying closure can accommodate 18 inde- 
pendent functions. Could this case, for instance, correspond to the emergence 
of two light cones, each with its 9 independent components? 

Unfortunately, a computation of the Fresnel tensor and the Fresnel equation 
using the general asymmetric solution (|297|) to (|300|) of the closure relation does 
not help much to recognize a possible double light cone structure. The Fresnel 
equation is still a quartic equation, in general. 

In what follows, we will study a particular case in order to try to get an 
idea of the qualitative properties that an additional skewon piece can induce. 
More details can be found in |61J . We will consider the case in which K = 0, 
which implies V = 0, C = 00 and A = —B^ 1 . Furthermore, we decompose the 
arbitrary matrix B into its symmetric and antisymmetric parts, 

Bab = b ab + e abc n c , with b ab := B {ab) , n c := e cab B [ab] . (332) 

Note that b ab contributes only to ^'X an d n c only to ^x- In this case, by 
substituting ([332]) into (fl82|) . we obtain: 



2 



w= q ° 



(detb + n 2 ) 



q i-2ql(q 2 -( q n) 2 ) + (f + ( q n) 2 y . (333) 



Here we used the abbreviations n 2 := b ab n a n b , q 2 :— b q a q b , qn := q a n a , and 
b ab is the (symmetric) matrix of the minors of b ab . 

In general, the right hand side of (|333[) is neither a square of a quadratic 
polynomial nor a product of two quadratic polynomials. In other words, neither 
a light cone nor a birefringence (double light cone) structure arises generically. 



'i.e., a medium with no 'magneto-electrical' properties. 
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Formally, one can rewrite (|333|) as 



W 
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(detfc + n 2 ) 



(qn) 



(qn) 



(334) 



Thus, the question of the reducibility of the Fres nel e quation translates into 
the algebraic problem of whether the square root \/—q 2 is a real linear polyno- 
mial in q a . There are three cases, depending on the rank of the 3x3 matrix 
bab, see [61]. 

(i) When b a t has rank 3, no factorization into light cones is possible (the 
roots a are complex), unless n a = 0. This latter condition implies ( 2 'x = 0, and 
the previous results are recovered. 

(ii) When b a b has rank 2, i.e., detb = 0, but at least one of the minors is 
nontrivial. Then, without loss of generality, one can assume: 



bu bi2 
h 2 b 2 2 




(335) 



In order to avoid complex solutions, we have to assume that the minor b 



(X < 0, so that 



W 



/z<73. In this case we find 



b 11 (n 1 ) 2 + 2b 12 n 1 n 2 + b 22 (n 2 ) 2 
q 2 + (qin 1 + q 2 n 2 + q 3 (n 3 + fj,))' 

ql + (qin 1 + q 2 n 2 + q 3 (n 3 - fi)) 2 ) 



(336) 



i.e. two different light cones. We can read off, up to conformal factors, the 
components of the two corresponding 'metric' tensors defining the light cones: 



ii 
9x 



9i 



( 1 

(n 1 ) 2 nhi 2 

n x n 2 (n 2 ) 2 

\ n 1 (n 3 + fj) n 2 n 3 







n 1 ^ 3 + ji) 

2 3 

n n 
(n 3 + ^) 2 J 



(1 





V n 1 (n 3 





(n 1 ) 2 

1 2 

n n 







2\2 





n 1 (n 3 - 



M) 



(n 2 ) 



(337) 



(338) 



{n 3 -n? J 



One verifies that det(g\ j ) = det(g%) = (n 1 ) 2 (n 2 ) 2 b 33 = -(n x ) 2 (n 2 ) 2 ^ 2 < 0, so 
that both metrics have the correct Lorentzian signature. 

(iii) When the 3x3 matrix b a b has rank 1. In this case corresponds to the 
case 2 with fi = 0. The Fresnel equation reduces to a single light cone, but the 
resulting metric is degenerated, since det(<7 y ) = 0. 
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Thus, we have demonstrated with a special case that some of the asymmetric 
solutions of the closure relation can yield birefringence. 

We saw that the conditions of closure and symmetry of x are sufficient for 
the existence of a well-defined light cone structure. If any of these conditions 
is violated, the light cone structure seems to be lost. The necessary conditions 
have still to be found. 

5 Conclusions and Prospects 

In this work, we have developed a general framework for describing classical 
electrodynamics in a general 4-dimensional medium. We have paid special at- 
tention to identify the structures which are metric-independent, and which can 
therefore be applied to a great variety of particular cases, from generalized mod- 
els describing the electromagnetic properties of spacctime to classical optics in 
material media. 

One of the central results of this work is summarized in the Fresnel equation 
(|181[) and in particular in the Fresnel tensor (| 180|) . This equation determines 
the local properties of the propagation of waves, in particular it describes the 
geometry of the wave covectors. It is a remarkable result that this important 
equation could be derived for any linear medium. Furthermore, the result can 
also be used to describe the effective properties of electromagnetic perturbations 
in nonlinear media. Its structure, i.e., its dependence on the constitutive tensor 
of the corresponding medium is highly nontrivial. Cubic structures are not very 
common in physics. To the best of my knowledge, no previous derivation of the 
Fresnel equation has been given which is as general as that in section HOI see 
also the early work of Tamm . 

It is also important to emphasize that the Fresnel equation and the whole 
formalism is generally covariant. This means that no artificial or particular 
coordinate choices are necessary and that only quantities describing intrinsic 
properties, in our case the constitutive tensor, of the physical system under 
consideration enter in the formalism. Of course, inside a given medium, for 
example, one can choose specific coordinates which are useful for concrete cal- 
culations. This is usually the case if the medium possesses some symmetry, as 
defined in section [2. 7. II Then an adapted coordinate system can be used which 
exploits this symmetry, so that calculations become simpler. This is, however, 
only a convenient choice, determined by the properties of the system (the con- 
stitutive tensor), but not a necessary a priori ingredient. In a material medium 
which is not isotropic, for instance, there is no a priori reason to use cartesian 
coordinates, the later are rather coordinates adapted to an isotropic medium, 
as for instance the minkowskian vacuum. Therefore, it is highly satisfactory to 
be able to describe the electromagnetic properties of a material medium in a 
generally covariant way. 

We also studied the properties of the three different irreducible pieces which 
a general constitutive tensor can contain. We have proved that the first, the 
symmetric piece ^x is essential for the medium to admit well behaved wave 
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properties. This piece is therefore the most important object determining the 
geometry of wave covectors and the light cone structure under some particular 
conditions. The second piece ^X was shown to describe dissipative properties 
of the medium. A nonvanishing piece can lead to dissipation of electro- 
magnetic energy, even is the material's properties are time independent. Fur- 
thermore, the skewon piece does influence the wave propagation and therefore 
the light cone structure. We have shown, studying some particular cases, that 
the skewon field can even give rise to a double light cone structure. The third 
possible irreducible piece ( 3 'x of a constitutive tensor, the abelian axion, is more 
elusive. It does not contribute to the energy-momentum density of the electro- 
magnetic field and does not affect the local properties of the electromagnetic 
waves and thus not the light cone structure either. It is possible, however, that 
it could influence the properties of wave propagation over finite (long) distances. 
Examples of all the three irreducible pieces of a constitutive can be found in the 
literature. The skewon field has been discussed so far only for material media 
violating CP symmetry, but up to now not as a model for the classical electro- 
magnetic properties of spacetime. If one follows the common assumption that a 
fundamental physical system should be describable in terms of an action prin- 
ciple, then the detection of a skewon piece could be interpreted as indicating an 
underlying substructure of the system under study. Of course, we do not claim 
here that the skewon field is actually nonvanishing in spacetime. We do believe, 
however, that it is a possible piece which can be used to quantify well defined, 
in general dissipative, properties of spacetime. Additionally, the skewon piece 
can be used to model possible effects violating local Lorentz invariance. Within 
the framework of GR, where a metric is present, one can show that the skewon 
piece necessarily violates local Lorentz invariance, see appendix [Cl 

Furthermore, we have studied the conditions under which a light cone struc- 
ture is induced. We have seen that, if the constitutive tensor is symmetric and 
satisfies a closure relation, a light cone is, in fact, induced. We also saw that 
the necessary and sufficient conditions for a constitutive tensor to be written 
as proportional to the Hodge dual operator of some metric, are symmetry and 
closure. This one-to-one relationship is valid when wc formulate it in terms of 
the constitutive tensor. However, it is still an open question whether symmetry 
and closure are necessary conditions for defining a light cone, in the sense that 
the quartic Fresnel equation reduces to a quadratic equation for the wave covec- 
tors, as discussed at the beginning of section [U We proved that symmetry and 
closure are sufficient conditions for such a reduction. Strictly speaking, closure 
of the whole constitutive tensor is not a necessary condition since an additional 
axion piece does not disturb the Fresnel equation. However, it may be that 
closure of only the principal irreducible piece ^'X an( i symmetry, i.e. ( 2 )x = 0, 
could be the necessary conditions. A proof of this conjecture is, however, still 
missing. Non-trivial examples supporting this conjecture can be found in [55] , 

Other interesting open issues which could be investigated are: We have seen 
that the constitutive tensor defines a certain generalization of the conformal 
properties of spacetime. Does it also defines a (generalized) affine structure? To 
answer this question, it would be interesting to study some additional properties 
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of the propagation of waves, for instance, how the polarization vector propagates 
along a light ray. It should be possible to find an answer to this question within 
our general formalism and it is expected to involve some affine properties (under 
which parallel displacement is the polarization vector constant along the light 
ray?). 



A The situation so far 
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Figure 8: Different particular cases for the constitutive tensor and the cor- 
responding behavior of the Fresnel equation and, therefore, of the light cone 
structure . The quantities K and b refer to the 3-dimensional matrices in sec- 
tion ESD 



B Electrodynamics in a material medium 

In a material medium, treated as a continuum, one can describe the properties 
of electromagnetic fields by means of macroscopic Maxwell's equations. A ma- 
terial medium can be defined in general terms as a region with a given total 
charge at macroscopic scales, but with a microscopic charge substructure. This 
substructure corresponds to the individual electrons, ions, etc. which form the 
medium. The particular distribution and dynamics of the charges forming the 
material are in general unknown or difficult to model. Therefore, it is useful to 
describe the electrodynamics of materials only in term of the so called 'external' 
charges, which are those that can, in principle, be manipulated in experiments. 
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The procedure to define macroscopic Maxwell equations for the material 
medium starting from the microscopic ones consist then in the separation of 
the total charge in the sum of the two contributions originating 'from the inside' 
(bound or polarizational charge) and 'from the outside' (free or external charge): 

J = J mat + J 0Xt . (339) 

Bound charges and bound currents are inherent characteristics of matter deter- 
mined by the medium itself. In contrast, external charges and external currents 
appear in general outside and inside matter and can be prepared for a specific 
purpose. 

The external charge is assumed to be conserved, and consequently the bound 
charge is conserved, too: 

dJ mat = 0j dJ cxt = Q ( 34 q) 

As in the case of vacuum, see Refs. [52] [24] [27], this allows us to introduce the 
corresponding excitation i^ mat as a 'potential' for the bound current: 

dH mat = J mat . (341) 

In a (1+3) decomposition, the 6 components of H m£Lt are identified with the 
polarization P and magnetization M. 

Defining now the external excitation if cxt as 

ff oxt := H - ff mat , (342) 

we find the inhomogeneous Maxwell equation 

diJ cxt = J oxt . (343) 

Here, only external quantities are involved. 

It remains to specify the constitutive law. In addition to the spacetime re- 
lation, H = H(F), the knowledge of the internal structure of a medium yields 
the macroscopic excitation 7J mat (i.e., the polarization and magnetization) as a 
function of the electromagnetic field strength F (and possibly of some thermo- 
dynamical variables describing the material continuum) . Then the constitutive 
law of the material is given by 

iJ oxt = H cxt {F) = H{F) - iJ mat (F). (344) 

C Local Lorentz invariance 

Given a lorentzian metric, we can define the notion of local Lorentz invariance. 
Let T il '"' lp be the contravariant coordinate components of a tensor field and 
T Ql - Qp := e^" 1 ■ ■ ■ a a " T 11 —' 1 " its frame components with respect to an or- 
thonormal frame e a — e l a eV A tensor is said to be local Lorentz invariant at a 
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given point, if its frame components are invariant under a local Lorentz rotation 
of the orthonormal frame. Similar considerations extend to tensor densities. 

There are only two geometrical objects which are numerically invariant under 
local Lorentz transformations: the Minkowski metric r) a p and the Levi-Civita 
tensor density e a ^s- Thus 

T a ^ s = tf>(x) (?7 Q V 5 " r/% Q<5 ) + a(x) e afi ^ (345) 

is the most general form of the frame components of a locally Lorentz invariant 
contravariant fourth rank tensor with the symmetries T^ kl = —T^ lkl = _7"w«_ 
Here (f> and a are scalar and pseudo-scalar fields, respectively. Therefore, back 
in coordinate components, we find that 

X iiU = <t>{x) V=9 (ffV ~ 9 ]k 9 a ) + a(x) e^ kl (346) 

is the most general form of a locally Lorentz invariant constitutive tensor allowed 
in a space with a lorentzian metric. Notice that this constitutive tensor is 
necessarily symmetric. In other words, any additional asymmetric piece ~ ^X 
would violate local Lorentz invariance. This was found before in the particular 
examples of constitutive tensors studied by Nieves and Pal, see section 12.6.31 
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